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Control  of  vortex  breakdown  continues  to  be  of  vital  importance  because  the 
breakdown  may  have  a  considerable  effect  on  aircraft  performance,  such  as  the  effect  on  the 
time-averaged  lift  force  (Lee  and  Ho,  1990).  Breakdown  is  also  important  because  of  the 
unsteady  nature  of  flow  downstream  of  vortex  breakdown.  The  imsteadiness  may  affect  the 
stability  of  the  aircraft  and  also  cause  buffeting.  The  objectives  of  this  work  are  to 
understand  the  physics  of  tmsteady  flow  phenomena  in  vortex  breakdown  region,  and  to 
develop  techniques  for  active  (feedback)  and  passive  (open  loop)  control  of  vortex 
breakdown  over  steady  and  unsteady  delta  wings. 

The  control  problem  is  particularly  important  for  unsteady  delta  wings  as  well  as  for 
unsteady  free  stream  (LeMay  et  al.,  1990;  Gursul  and  Ho,  1994).  The  vortex  breakdown 
may  appear  on  the  wing  even  if  no  breakdown  is  observed  in  the  steady  case  (Gursul  and 
Ho,  1993).  The  time-dependent  nature  of  breakdown  location  may  affect  imsteady  loading 
on  the  wing.  One  of  the  objectives  of  this  work  is  to  study  the  possibility  of  controlling  the 
breakdown  in  imsteady  flows. 

It  is  well  known  that  there  are  two  important  parameters  which  determine  the 
breakdown  location:  swirl  angle  and  external  pressure  gradient  outside  the  vortex  core. 
The  principle  of  any  control  method  is  to  alter  either  or  both  of  these  parameters.  Blowing 
and  suction  in  the  tangential  direction  along  a  rounded  leading-edge  (Wood  et  al.,  1990;  Gu 
et  al.,  1993)  and  suction  applied  around  the  vortex  axis  (Werle,  1960;  Parmenter  and 
Rockwell,  1990)  were  shown  to  delay  breakdown.  The  former  is  believed  to  be  due  to  a 
change  in  the  strength  and  location  of  the  leading-edge  vortices.  This  is  achieved  by 
affecting  the  location  of  separation  on  the  rounded  leading-edge.  Direct  control  of  the 
separation  point  modifies  the  vortex  properties  such  that  a  reduction  in  the  swirl  angle  is 
achieved.  The  purpose  of  applying  suction  around  the  vortex  axis  is  to  reduce  the  local 
adverse  pressure  gradient,  which  is  achieved  by  accelerating  the  axial  flow  along  the  core. 

Since  all  of  the  vorticity  of  the  leading-edge  vortices  originates  from  the  separation 
point  along  the  leading-edge,  leading-edge  devices  are  particularly  attractive  tools  that  can 
be  used  to  influence  the  strength  and  structure  of  these  vortices.  For  example,  leading-edge 
flaps  are  known  to  be  capable  of  controlling  the  circulation  and  location  of  the  leading-edge 


3 


vortices  (Spedding  et  al.,  1987;  Karagounis  et  al.,  1989;  Cheng  et  al,  1988).  These  studies 
concentrated  primarily  on  the  effect  of  leading-edge  flaps  at  low  angle  of  attack  where 
vortex  breakdown  was  not  observed  over  the  wing.  Recently,  the  effect  of  vortex  cavity 
flaps  at  high  angles  of  attack  has  been  reported  (Schaeller  et  al.,  1994).  One  of  the  purposes 
of  this  work  is  to  explore  the  effect  of  leading-edge  flaps  on  vortex  breakdown  phenomena. 
The  stationary  as  well  as  oscillating  flaps  may  provide  further  insight  into  the  possibility  of 
control  of  vortex  breakdown.  We  also  studied  vortex  breakdown  over  a  delta  wing  with 
variable  sweep  angle  since  this  method  allows  a  direct  control  of  circulation  of  leading-edge 
vortex.  Again,  periodic  variations  of  sweep  angle  was  studied  with  the  control  of  vortex 
breakdown  in  mind.  Application  of  this  technique  in  unsteady  flows  was  demonstrated  for 
a  pitching  delta  wing. 


2.  EXPERIMENTAL  FACILITIES 

Experiments  were  carried  out  in  a  wind  tunnel  and  a  water  channel,  both  with  a 
cross-sectional  area  of  61  cm  by  61  cm.  The  first  model  delta  wing  (see  Figure  1)  was  only 
used  in  wind  tunnel  testing.  It  had  a  sweep  angle  of  70°.  It  was  made  of  3/4in.(19mm) 
thick  Plexiglas  and  had  a  chord  length  of  c=14  in.(356mm).  The  surface  was  flat  whereas 
the  leading  edges  were  beveled  at  45°  on  the  windward  side.  On  the  lee  surface  grooves 
were  made  for  purposes  of  pressure  measurements.  The  grooves  were  covered  by 
aluminum  plates  on  which  pressure  taps  of  1  mm  diameter  and  a  depth  of  2  mm  were 
located.  The  plates  were  flush  mounted  to  the  surface  so  that  no  interference  of  flow  would 
be  introduced.  Unsteady  pressure  measurement  locations  in  the  streamwise  direction  were 
along  the  centerline  and  along  the  ray  y/s=0.5  from  the  apex.  Additional  pressure  taps  for 
mean  pressure  measurements  were  located  in  the  spanwise  direction  at  x/c=0.6  and  x/c=08. 
All  unused  taps  were  sealed  by  tape  during  the  experiment.  The  Reynolds  number  was 
varied  from  127,000  to  254,000,  although  most  of  the  results  are  for  Re=254,000.  At  the 
largest  angle  of  attack,  the  blockage  ratio  was  0.08. 

The  other  two  model  delta  wings,  one  with  variable  sweep  angle  and  one  with 
leading-edge  flaps,  are  shown  in  Figures  2  and  3,  respectively.  For  both  models,  the  base 


model  has  a  sweep  angle  of  A=70°.  The  chord  lengths  are  c=254  mm  and  268  mm  for  the 
model  with  flaps  and  the  variable  sweep  delta  wing,  respectively.  The  Reynolds  number 
was  around  Re=50,000  for  the  water  channel  experiments.  The  flap  angle  5  (measured  from 
the  upper  surface)  could  be  varied  up  to  180°.  For  the  variable  sweep  delta  wing,  the  range 
of  sweep  angle  was  from  A=60°  to  A=70°.  The  base  plate  of  the  models  was  made  of 
Plexiglas  and  the  thdn  plates  were  made  of  stainless  steel.  The  motion  of  the  flaps  was 
guided  with  gears  and  joints  which  were  connected  to  a  ring.  The  flaps  were  opened  by 
rotating  this  ring.  The  rotation  of  the  ring  was  achieved  by  pulling  a  cable  which  was 
attached  to  the  ring.  For  the  variable  sweep  delta  wing,  the  two  plates  were  used  to  change 
the  sweep  angle.  The  motion  of  the  plates  was  guided  with  a  cable-pulley-pin  system.  The 
wing  was  opened  by  pulling  the  cable.  For  both  models,  a  spring  assured  the  closing  of  the 
wing  or  flaps  when  the  tension  on  the  cable  was  released.  A  bicycle  brake  cable  whose 
sheath  was  fixed  near  the  trailing-edge  and  outside  the  water  channel  was  attached  to  a 
drum  which  was  driven  by  a  DC  motor  position  control  system  (see  Figure  4).  A  sine  wave 
from  a  function  generator  was  used  to  drive  the  DC  motor  position  control  system.  A 
potentiometer  was  used  to  monitor  the  variation  of  the  flap  angle  or  sweep  angle.  This 
flexible  system  allowed  the  use  of  flaps  or  variable  sweep  even  for  a  pitching  motion  of  the 
delta  wings. 

The  experiments  were  conducted  both  for  stationary  and  pitching  delta  wings.  The 
pitching  mechanism  was  similar  to  the  one  used  by  LeMay  et  al  (1990).  A  displacement 
transducer  was  used  to  monitor  the  variation  of  the  angle  of  attack.  A  variable  speed  DC 
motor  and  a  speed  controller  were  used  to  drive  the  pitching  mechanism.  For  the  combined 
motion  of  pitching  and  harmonic  variations  of  sweep  angle,  the  input  signal  to  the  DC 
motor  position  control  system  was  obtained  by  applying  a  phase  shift  to  the  displacement 
transducer  signal  from  the  pitching  mechanism.  This  provided  the  desired  phase  angle 
between  angle  of  attack  a(t)  and  sweep  angle  A(t). 

Flow  visualization  of  vortex  breakdown  was  done  by  injecting  fluid  with  food 
coloring  dye  near  the  apex  of  the  models.  The  flow  visualization  was  videotaped  for  further 
analysis.  The  velocity  was  measured  with  a  single  component  laser  Doppler  velocimetry 


(LDV)  system.  The  measurement  uncertainties  for  breakdown  location  (x/c)  and  mean 
velocity  were  estimated  as  0.4%  and  1%,  respectively.  Ensemble-averaging  technique  was 
applied  to  the  velocity  signals  in  the  case  of  periodic  variations  of  flap  angle  or  sweep 
angle.  Typically,  100  cycles  were  averaged  in  order  to  obtain  the  phase-averaged  velocity. 
The  standard  deviation  from  the  phase-averaged  value  was  also  calculated. 

3.  RESULTS 

3.1.  Unsteady  Flow  Phenomena  in  Vortex  Breakdown  Region 

In  the  wind  tunnel  experiments  with  the  first  model  delta  wing,  coherent  pressure 
fluctuations  have  been  detected  on  the  upper  wing  surface  as  the  angle  of  attack  increases 
and  breakdown  location  moves  over  the  wing  as  shown  in  Figure  5.  The  time  traces  shown 
in  the  figure  were  taken  at  a  location  close  to  trailing-edge  (x/c=0.93)  and  approximately 
underneath  the  leading-edge  vortex  (y/s=0.5).  At  the  smallest  angle  of  attack  a=l5°,  there 
is  no  breakdown  over  the  wing.  At  larger  angles  of  attack,  the  breakdown  moves  over  the 
wing  causing  increased  pressure  fluctuations. 

Simultaneous  pressure  measurements  at  two  locations  in  the  spanwise  direction 
indicated  that  the  disturbances  rotate  in  the  same  direction  as  the  leading-edge  vortex  (see 
Figure  6).  Moreover,  simultaneous  hot-wire  measurements  at  two  points  which  are  at  either 
side  of  the  vortex  axis  showed  that  there  existed  a  phase  angle  of  approximately  180  deg. 
between  the  signals  (Figure  7).  This  implies  that  the  disturbances  are  helical  waves  with 
azimuthal  wavenumber  n=l,  if  the  instability  is  presented  as  exp{i(<xx+n^-(i)t)}.  It  was 
shown  that  this  type  of  unsteady  flow  is  a  result  of  the  instability  of  the  time-averaged 
swirling  flow  with  a  wake-like  axial  profile  (Garg  and  Leibovich,  1979;  Gursul,  1994),  as  is 
found  in  breakdown  regions,  and  is  observed  in  other  swirling  flows  as  well. 

Measurements  showed  that  the  dominant  frequency  detected  from  the  pressure 
fluctuations  on  the  upper  surface  is  not  constant,  but  vary  with  the  streamwise  distance. 
Other  main  parameters  of  the  instability  such  as  the  wavelength  and  phase  speed  were 
determined  for  a  large  range  of  angle  of  attack  (Gursul,  1994,  see  Appendix  A). 


3.2.  Time  Scales  of  Vortex  Breakdown 

It  has  been  observed  in  several  experiments  that  vortex  breakdown  location  is  not 
steady  and  exhibits  fluctuations  in  the  stream  wise  direction  (Lowson,  1964;  Garg  and 
Leibovich,  1979;  Payne  et  al.,  1988;  Thompson  et  al.,  1991).  One  possible  source  is  the 
hydrodynamic  instability  of  the  “breakdown  wake”  (i.e.,  the  flow  in  the  wake  of  the  vortex 
breakdown)  where  wake-like  axial  velocity  profiles  are  observed.  The  quasi-periodic  nature 
of  the  flow  in  the  wake  of  breakdown  may  interact  with  the  breakdown  process  through 
upstream  influence.  It  has  not  been  possible  to  determine,  from  the  previous  experiments, 
the  role  of  the- helical  mode  instability  in  the  imsteady  nature  of  breakdovra  location. 

We  carried  out  flow  visualization  and  laser  Doppler  velocimetry  (LDV) 
measurements  in  the  water  channel.  The  results  (see  Figure  8)  indicated  that  the 
fluctuations  of  vortex  breakdown  location  occur  at  much  lower  frequencies  than  the 
frequency  of  the  hydrodynamic  instability  of  the  flow  in  the  wake  of  breakdown.  It  is  thus 
suggested  that  the  helical  mode  instability  does  not  influence  the  unsteady  nature  of 
breakdown  location  (Gursul  and  Yang,  1995,  see  Appendix  B). 

In  order  to  understand  the  nature  of  the  organized  fluctuations,  simultaneous  flow 
visualization  of  the  two  leading  edge  vortices  was  performed  for  A=70°  and  a=37°.  The 
cross-spectrum  of  the  fluctuations  for  the  two  breakdovras  (the  left  and  right  vortices)  is 
shown  in  Figure  9.  It  is  obvious  that  both  breakdowns  oscillate  at  the  same  frequency  and 
the  phase  angle  between  them  is  approximately  180°.  In  other  words,  there  is  an 
antisymmetric  motion  of  breakdown  locations  for  left  and  right  vortices.  There  is  clearly  an 
interaction  between  the  opposite  leading-edge  vortices.  However,  the  exact  nature  of  this 
interaction  is  not  known. 

As  a  summary,  the  spectrum  of  unsteady  flow  phenomena  over  delta  wings  as  a 
function  of  dimensionless  frequency  is  shown  in  Figure  10.  In  this  graph,  the  range  of  the 
dominant  frequency  of  the  oscillations  of  breakdo-wn  location  was  obtained  from  the  flow 
visualization  for  different  angles  of  attack.  The  range  of  the  helical  mode  instability  was 
available  from  the  previous  experiments  (Gursul,  1994).  The  frequency  range  for  the 
Kelvin-Helmholtz  instability  of  the  shear  layer  separated  from  the  leading-edge  was 


estimated  from  the  experiments  and  numerical  simulations  (Gad-el-Hak  and  Blackwelder, 
1987;  Gordnier  and  Visbal,  1994).  When  compared  with  the  frequency  of  the  helical  mode 
instability,  the  frequency  range  of  the  spectrum  of  breakdown  location  for  a  stationary  delta 
wing  is  much  closer  to  the  frequency  range  of  typical  aerodynamic  maneuvers  (up  to 
fc/Uoo«0.03).  The  response  of  breakdown  location  and  possible  coupling  between  the  wing 
motion  and  breakdown  location  in  this  frequency  range  is  very  important. 

3.3.  Active  (Feedback)  Control  of  Vortex  Breakdown 

The  pressure  fluctuations  induced  by  the  helical  mode  instability  of  vortex 
breakdown  can  be  used  to  control  the  vortex  breakdown  location.  As  shown  in  Figure  5, 
with  the  increasing  angle  of  attack  as  the  breakdown  moves  toward  the  apex,  the  amplitude 
of  pressure  fluctuations  increases.  Hence  the  rms  value  of  pressure  was  chosen  as  the 
control  variable.  The  basic  idea  of  the  breakdown  control  is  to  alter  the  circulation  of  the 
leading-edge  vortex  as  the  control  variable  varies.  This  has  been  achieved  with  the  model 
delta  wing  with  variable  sweep.  A  feedback  control  loop  was  used  to  control  the  amplitude 
of  pressure  fluctuations  at  a  single  point  on  a  delta  wing.  The  block  diagram  of  the 
feedback  control  system  is  shown  in  Figure  11.  Further  information  is  available  in  the 
published  article  (Gursul  et  al.,  1995,  see  Appendix  C). 

First,  the  active  control  experiments  were  conducted  for  a  stationary  delta  wing  at 
angle  of  attack  a=23  deg.  (see  Figure  12).  Initially  the  sweep  angle  was  set  to  A=60  deg. 
For  this  value  of  sweep  angle,  the  vortex  breakdown  location  is  over  the  wing  and  closer  to 
the  apex.  Once  the  control  loop  is  activated,  it  varies  the  sweep  angle  dynamically  until  the 
amplitude  of  pressure  fluctuations  are  below  a  threshold  value  (or  until  the  breakdown 
location  moves  downstream  of  trailing  edge).  In  figure  12,  the  plots  of  pressure  fluctuation, 
feedback  voltage,  and  sweep  angle  are  shovra.  Around  t=5  seconds,  the  feedback  control 
was  turned  on.  It  is  seen  that  the  sweep  angle  increases  and  becomes  nearly  constant 
around 
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A=68.5  deg.  after  a  small  overshoot.  This  new  sweep  angle  corresponds  to  a  case  where  the 
breakdown  location  is  downstream  of  the  wing.  The  decrease  in  pressure  fluctuations  and 
feedback  voltage  accompany  the  variation  of  sweep  angle. 

The  active  control  experiments  for  the  pitching  delta  wing  were  conducted  in  a 
similar  way.  As  the  delta  wing  was  pitched  periodically  between  15  deg  and  23  deg,  the 
initial  sweep  angle  was  set  to  60  deg.  For  this  range  of  angle  of  attack,  the  vortex 
breakdown  is  over  the  wing  during  the  whole  cycle.  The  plots  of  angle  of  attack,  pressure 
fluctuations,  feedback  voltage  and  sweep  angle  are  shown  in  Figure  13.  When  the  control 
loop  is  turned  on  suddenly,  the  sweep  angle  increases  and  becomes  nearly  constant  around  a 
new  value  for  which  the  breakdown  location  is  downstream  of  the  wing. 

6.4.  Vortex  Breakdown  over  a  Pitching  Delta  Wing. 

Several  experimental  studies  of  vortex  breakdown  over  pitching  delta  wings  have 
shown  that  there  is  a  phase  shift  between  the  motion  of  the  wing  and  the  movement  of 
breakdown  location  (LeMay  et  al.,  1990;  Wolffelt,  1987;  Woodgate,  1971;  Atta  and 
Rockwell,  1990).  The  phase  delay  mainly  depends  on  the  reduced  frequency  k=(i)c/2U„ 
(LeMay  et  al.  1990)  and  can  be  as  large  as  180°.  Nevertheless,  the  physics  of  this  phase  lag 
is  not  well  understood. 

It  has  been  suggested  that  delta  wing  vortex  breakdown  in  unsteady  free  stream  is 
driven  by  the  external  pressure  gradient  (Gursul  and  Ho,  1994).  Likewise,  the  observed 
phase  delay  of  the  breakdown  location  for  pitching  delta  wing  may  be  related  to  the 
variations  of  the  adverse  pressure  gradient  on  the  wing  surface.  In  order  to  test  this 
hypothesis,  unsteady  pressure  on  a  pitching  delta  wing  was  measured  and  correlated  with 
the  movement  of  the  breakdown  location. 

It  was  shown  earlier  that  vortex  breakdown  flowfields  over  delta  wings  show 
helical-mode  instability  which  produces  coherent  pressure  fluctuations  on  the  suction 
surface.  With  increasing  angle  of  attack,  the  amplitude  of  the  pressure  fluctuations  due  to 
the  signature  of  breakdown  increases.  When  the  delta  wing  is  pitched  sinusoidally 
(between  29°  and  39°),  the  instantaneous  pressure  signature  of  the  breakdown  varies  as  the 


breakdown  location  moves  back  and  forth  in  the  streamwise  direction  (Figure  14).  It  is 
evident  that  the  maximum  fluctuation  level  due  to  the  breakdown  occurs  not  at  the 
maximum  angle  of  attack,  but  later,  at  a  smaller  angle  of  attack  as  the  wing  pitches  down. 
In  order  to  quantify  the  phase  delay,  the  phase-averaging  technique  was  applied  to  the 
pressure  signal,  i.e. 

Cp  =  (Cp}  +c/ 

where  the  <Cp>  is  the  phase-averaged  pressure  coefficient.  The  standard  deviation  from  the 
phase-averaged  value. 


was  calculated  as  a  function  of  angle  of  attack.  The  variation  of  the  standard  deviation 
shown  in  Figure  15  reveals  hysteresis  loops.  The  loops  become  wider  with  increasing 
reduced  frequency,  which  is  an  indication  of  increasing  phase  delays. 

In  order  to  investigate  the  possible  relationship  between  the  pressure  field  and  the 
dynamic  character  of  breakdown  location,  pressure  measurements  were  carried  out  along 
the  centerline  (y/s=0)  and  y/s=0.5.  For  the  leading-edge  vortex,  the  pressure  gradient  on  the 
suction  surface  of  the  wing  acts  as  the  external  pressure  gradient.  However,  interpretation 
of  the  pressure  along  the  y/s=0.5  is  not  straightforward,  since  the  measurement  locations  are 
approximately  underneath  the  vortex.  Whether  the  unsteady  variation  of  pressure  is  the  end 
result  of  the  vortex  breakdown  remains  uncertain.  On  the  other  hand,  along  the  centerline 
y/s=0,  the  effects  of  vortex  on  pressure  are  negligible.  This  is  verified  by  a  simple  test.  A 
control  cylinder  with  diameter  d/c=0.047  was  placed  at  x/c=0.92,  z/c=0.37  (for  a=29°) 
which  caused  a  premature  breakdown.  The  mean  pressure  coefficient  along  the  spanwise 
direction  at  x/c=0.80  for  the  cases  with  and  without  the  cylinder  shows  that  there  is  a  large 
reduction  in  the  induced  pressure  underneath  the  vortex  due  to  the  breakdown,  whereas  the 
effects  are  minimum  at  the  centerline  (Figure  16). 

For  the  pitching  motion,  the  variation  of  the  phase-averaged  pressure  along  the 
centerline  as  a  function  of  time  is  shown  in  Figure  17  for  two  different  values  of  reduced 
frequency.  The  angle  of  attack  varies  as 

a  =  34°  -  5°  cos  ©t. 


10 


The  pressure  along  the  centerline  varies  in  magnitude  as  a  function  of  time.  For  the  smaller 
reduced  frequency  (k=0.19),  the  maximum  pressure  (suction)  occurs  when  the  angle  of 
attack  is  closer  to  the  maximum  angle  of  attack;  however,  there  exists  a  time  delay.  For 
k=0.35,  this  phase  delay  is  larger.  The  pressure  field  for  a  pitching  wing  seems  to  be 
delayed  in  time  compared  to  the  quasi-steady  case.  This  phase  delay  is  close  to  that  of  the 
breakdown  location  obtained  from  pressure  fluctuations  (i.e.,  from  the  variation  of  the 
standard  deviation  shown  in  Figure  15).  Thus,  there  exists  a  definite  relationship  between 
the  pressure  field  on  the  wing  and  the  location  of  the  vortex  breakdown. 

The  Fourier  analysis  of  the  phase  averaged  signal  showed  that  most  of  the  energy 
was  concentrated  at  the  fundamental  frequency.  The  phase  angle  between  the  angle  of 
attack  a  and  (-Cp)  was  calculated  at  each  measurement  station  and  plotted  in  Figure  18. 
The  phase  delay  increases  with  increasing  reduced  frequency.  At  large  reduced  frequencies, 
there  is  also  considerable  variation  of  the  phase  angle  in  the  streamwise  direction. 
However,  it  should  be  kept  in  mind  that  the  amplitude  of  the  pressure  variations  becomes 
smaller  toward  the  trailing-edge  where  the  largest  phase  angles  are  observed. 

The  variation  of  phase  lag  between  the  wing  motion  and  <-Cp>  is  shown  as  a 
function  of  reduced  frequency  in  Figures  19.  Also  shown  are  the  phase  lags  obtained  from 
the  pressure  fluctuations  in  addition  to  those  determined  by  flow  visualization  for  a  similar 
model.  It  is  evident  that  there  exists  a  definite  relationship  between  the  pressure  field  on  the 
wing  surface  and  the  location  of  the  vortex  breakdown. 

In  conclusion,  variation  of  the  phase-averaged  pressure  along  the  wing  centerline  as 
a  function  of  time  shows  that  there  is  a  phase  delay  between  the  wing  motion  and  the 
pressure  field.  This  phase  delay  is  close  to  the  phase  delay  of  the  movement  of  breakdown 
location.  The  phase  delay  of  pressure  increases  with  increasing  reduced  frequency  and  is 
not  sensitive  to  Reynolds  number  variation.  Similar  phase  lags  of  pressure  were  found  for 
smaller  angle-of-attack  range  for  which  the  breakdown  is  absent.  This  confirms  that  the 
observed  time  delays  in  breakdown  location  is  related  to  the  variation  of  the  external 
pressure  gradient  generated  by  the  wing.  Further  information  can  be  found  in  the  published 
journal  article  (Gursul  and  Yang,  1995,  see  Appendix  D). 
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3.5.  Control  of  Vortex  Breakdown  with  Leading-edge  Flaps 

The  location  of  vortex  breakdown  as  a  function  of  the  flap  angle  5  is  shown  in 
Figure  20  for  several  values  of  angle  of  attack.  At  low  angles  of  attack  (a=16°  and  a=20°), 
as  the  flap  angle  decreases  from  1 80°,  the  breakdown  location  moves  toward  the  trailing- 
edge.  This  effect  is  similar  to  that  of  increasing  sweep  angle  (or  decreasing  aspect  ratio). 
The  breakdown  location  suddenly  moves  in  far  wake  as  it  approaches  the  trailing-edge.  At 
larger  angles  of  attack  (a=25°  and  a=30°),  the  breakdown  location  does  not  change  much 
initially  as  the  flap  angle  decreases  from  180°,  but  it  moves  rapidly  toward  the  trailing-edge. 
At  the  largest  angle  of  attack  (a=35°),  the  breakdown  location  moves  downstream  first, 
with  decreasing  flap  angle,  and  then  upstream  for  5<65°. 

Obviously,  the  effect  of  flaps  on  vortex  breakdown  strongly  depends  on  the  angle 
of  attack.  If  one  considers  the  response  of  vortex  breakdown,  for  example  X5d/c>0.40,  there 
is  a  clear  difference  between  the  variations  of  breakdown  location  with  flap  angle  at  a=16°, 
20°  and  a=25°,  30°.  The  slopes  of  the  variations  are  very  different.  In  order  to  gain  further 
insight,  detailed  flow  structures  were  studied  for  a=20°  and  25°. 

A  variety  of  flow  patterns  were  observed  depending  on  flap  angle  and  angle  of 
attack.  As  an  example,  constant  levels  of  the  mean  axial  velocity  in  a  cross  plane  upstream 
of  breakdown  (x/c=0.40)  are  shovra  for  two  flap  angles  and  for  no  flaps  in  Figure  21.  The 
presence  of  the  flaps  produces  an  adverse  pressure  gradient  in  the  spanwise  direction  and 
increases  the  size  of  the  separated  flow  region  on  the  upper  surface  of  the  wing.  This  also 
causes  a  marked  change  in  the  structure  of  the  leading-edge  vortices.  For  6=60°,  the  core 
of  the  leading-edge  vortex  is  much  larger  compared  to  that  of  the  wing  with  no  flaps.  For 
6=100°,  much  larger  axial  velocities  are  observed  in  the  core.  This  also  implies  larger 
pressure  drop  at  the  core,  which  in  turn  is  a  result  of  the  larger  swirl  velocities.  Outside  of 
the  flaps,  either  high  speed  regions  or  flow  reversal  (due  to  the  separated  flow)  are  observed 
depending  on  flap  angle  and  angle  of  attack. 

Variation  of  swirl  velocity  was  also  investigated  by  traversing  across  the  vortex 
core.  The  maximum  swirl  velocity  is  shown  in  Figure  22  as  a  function  of  the  flap  angle. 
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The  arrows  on  the  figure  indicate  the  flap  angle  at  which  breakdown  is  at  the  trailing-edge 
of  the  wing.  It  seems  that  the  vortex  breakdown  appears  on  the  wing  when  the  swirl  level 
reaches  a  certain  value.  For  a=20°,  the  swirl  level  does  not  change  much  with  the  flap 
angle  6  when  the  breakdown  is  over  the  wing.  On  the  other  hand,  for  a=25°,  the  swirl  level 
increases  rapidly  with  the  flap  angle  (when  the  breakdown  is  over  the  wing),  explaining  the 
sensitivity  of  breakdown  location  -with  respect  to  the  variations  in  the  flap  angle. 

In  addition  to  the  changes  in  the  swirl  level  (or  strength  of  the  vortex),  the  location 
of  the  vortex  core  (in  the  cross  plane)  exhibits  variations  with  the  flap  angle.  The  location 
of  the  vortex  core  was  obtained  from  the  velocity  measurements  (not  shown  here),  which 
suggested  that  the  existence  of  the  flap  may  affect  the  pressure  gradient  outside  of  the 
vortex  core.  For  example,  for  a=20°,  there  is  a  gradual  variation  of  breakdown  location 
with  the  flap  angle,  although  the  swirl  level  does  not  change  much  in  the  same  range.  This 
implies  that  the  external  pressure  gradient  for  the  vortex  core  is  affected  by  the  changes  in 
the  geometry  of  the  wing. 

Experiments  with  oscillating  flaps  were  also  conducted.  The  variation  of  vortex 
breakdown  location  for  harmonic  variations  of  flap  angle  between  60°  and  180°  is  shown  in 
Figure  23  together  with  that  of  the  steady  flap  angles.  The  reduced  frequency  is 
k=coc/2U«,=0.4,  where  (o  is  the  radial  frequency  and  is  the  free  stream  velocity.  The 
variation  of  breakdown  location  reveals  a  hysteresis  loop.  This  is  an  indication  of  a  time 
lag  response,  which  is  similar  to  the  response  of  a  pitching  delta  wing. 

In  conclusion,  it  was  shown  that  the  effect  of  flaps  and  sensitivity  of  breakdown 
location  strongly  depends  on  angle  of  attack.  Velocity  measurements  showed  that,  with  the 
varying  flap  angle,  large  changes  take  place  in  the  vortex  core  diameter  and  location,  in 
axial  and  swirl  velocity  profiles,  and  in  the  secondary  vortex  structure.  In  general,  with  the 
increasing  flap  angle,  larger  swirl  velocities  are  observed.  The  variation  of  maximum  swirl 
velocity  seems  to  be  correlated  with  the  variations  of  breakdown  location  as  the  flap  angle 
is  varied.  In  addition,  experiments  with  oscillating  flaps  showed  that  there  exists  a 
hysteresis  loop  in  the  variation  of  breakdown  location.  This  indicated  a  time-lag  response 
similar  to  that  of  a  pitching  delta  wing. 


3.6.  Control  of  Vortex  Breakdown  with  Variable  Sweep 

In  static  experiments,  the  location  of  vortex  breakdown  moves  toward  the  apex  with 
decreasing  sweep  angle  (or  aspect  ratio),  as  known  from  a  number  of  experiments  with 
different  aspect  ratio  wings.  Therefore,  in  the  steady  case,  the  location  of  breakdown  shows 
a  monotonic  variation  with  the  sweep  angle,  as  opposed  to  the  variation  with  the  flap  angle. 
The  variation  of  breakdown  location  was  studied  for  harmonic  variations  of  sweep  angle, 
i.e. 

A=Ao+Ai  sin  cot 

The  variations  of  the  ensemble-averaged  breakdown  location  are  shown  in  Figure  24  for 
Ao=65°  and  Ai=2.5°  for  different  values  of  the  reduced  frequency,  which  show  hysteresis 
loops.  The  loops  become  narrower  with  increasing  reduced  frequency.  Actually,  at  very 
large  reduced  frequency  (k=2.0  and  k=4.0  (not  shown  here)),  the  ensemble-averaged 
breakdown  location  is  almost  constant.  The  phase  lag  of  vortex  breakdown  location  with 
respect  to  that  of  the  quasi-steady  case  was  found  by  plotting  the  time  history  of  breakdown 
location.  This  is  shown  in  Figure  25.  It  seems  that  the  phase  lag  reaches  a  maximum  at 
k=0.4  and  then  decreases.  This  time-lag  behavior  seems  to  be  an  inherent  response  of 
breakdown  location  in  unsteady  flows,  regardless  of  the  type  of  motion. 

In  order  to  understand  the  phase  lag  of  vortex  breakdown,  phase-averaged  velocity 
measurements  were  carried  out  in  a  cross  plane  (at  x/c=0.5)  upstream  of  breakdown 
location  for  a  reduced  frequency  of  k=0.4.  It  should  be  noted  that  the  flow  field  of  the 
leading-edge  vortex  develops  very  gradually  before  vortex  breakdown.  Indeed,  the  velocity 
changes  drastically  only  in  a  small  neighborhood  of  the  breakdown  location.  It  is  also 
known  that  the  swirl  level  of  the  leading  edge  vortex  does  not  vary  much  in  the  streamwise 
direction.  Therefore,  these  measurements  give  a  good  indication  of  the  upstream  conditions 
before  breakdown.  Contours  of  constant  phase-averaged  axial  velocity  are  shown  at 
different  sweep  angles  for  Ao=66°  and  A,=3°  in  Figure  26.  Also  shown  are  the  mean 
velocity  contours  in  the  steady  case  for  the  minimum  sweep  angle  A=63°  and  the  maximum 
sweep  angle  A=69°.  As  the  sweep  angle  varies  dynamically,  the  location  of  the  vortex  core 


(defined  as  the  location  at  which  the  phase-averaged  axial  velocity  is  maximum)  varies 
mostly  in  the  spanwise  direction.  Note  that  the  axial  velocity  contours  are  very  similar  for 
the  dynamic  and  static  cases.  Also  the  two  flow  fields  for  the  mean  sweep  angle  Ao=66° 
(for  increasing  and  decreasing  sweep  angle)  are  not  very  different  from  each  other.  All 
these  suggest  that  the  structure  of  leading  edge  vortex  is  not  affected  much  by  the 
unsteadiness  and  is  approximately  quasi-steady.  Nevertheless,  there  is  a  small  phase  lag  in 
the  development  of  the  vortical  flow.  This  is  evident  from  the  contours  of  constant  standard 
deviation  from  the  phase-averaged  velocity  shown  in  Figure  27.  It  should  be  noted  that  the 
shape  of  the  contours  of  the  standard  deviation  (or  rms  velocity  fluctuations  in  the  steady 
case)  is  related  to  the  main  vortex  as  well  as  the  shear  layer  and  secondary  separation 
region. 

The  variations  of  the  main  parameters  are  summarized  in  Figure  28.  The  variation 
of  sweep  angle  and  ensemble-averaged  breakdown  location  with  time  are  shown  at  the  top 
of  Figure  28  for  one  cycle  of  the  motion.  In  the  middle  of  the  figure,  the  location  of  the 
vortex  core  in  the  dynamic  and  static  cases  are  shown.  Note  that  the  phase  lag  of  the  core 
location  (with  respect  to  the  quasi-steady  case)  is  small  compared  to  the  larger  phase  lag  of 
breakdown  location  with  respect  to  the  sweep  angle  variations.  In  order  to  get  a  good 
indication  of  the  strength  of  the  vortex,  the  circulation  around  a  fixed  contour  (see  inset  in 
Figure  28)  was  calculated  from  the  phase-averaged  velocity  measurements.  Although  this 
contour  covers  the  main  vortex  and  most  of  the  secondary  vortex  and  shear  layer,  it  is  a 
reasonable  indicator  for  the  strength  of  the  main  vortex.  Moreover,  if  it  is  assumed  that 
conical  flow  exists,  the  circulation  increases  linearly  with  the  streamwise  distance  from  the 
apex.  In  this  case,  the  measured  circulation  at  a  fixed  streamwise  distance  is  related  to  the 
rate  of  increase  of  the  circulation  along  the  streamwise  direction  (or  the  rate  at  which 
vorticity  is  fed  into  the  leading-edge  vortex).  Note  that  the  phase  lag  of  the  circulation 
(with  respect  to  the  quasi-steady  case)  is  also  small  (bottom  of  Figure  28).  In  conclusion, 
the  unsteady  effects  on  the  vortex  core  location  and  strength  are  not  very  strong.  Therefore, 
the  large  phase  lag  of  breakdown  location  cannot  be  explained  with  the  variations  in  the 
upstream  conditions  since  the  location  and  strength  of  the  vortex  exhibit  smaller  phase  lags. 


This  suggests  that  the  variations  of  the  external  pressure  gradient  (which  is  one  of  the  two 
control  parameters  for  breakdown,  the  swirl  level  being  the  other  one)  should  be  responsible 
for  the  large  phase  lag  of  vortex  breakdown  location. 

Perturbations  applied  at  the  leading-edge  by  varying  sweep  angle  dynamically 
around  a  mean  value  may  provide  an  effective  way  for  breakdown  control  in  unsteady 
flows.  In  order  to  demonstrate  this,  a  pitching  motion  of  the  delta  wing  was  considered,  i.e. 
a=ao+ai  sin  cot 

The  well  known  hysteresis  loop  for  the  location  of  the  vortex  breakdown  is  shown  in 
Figure  29  for  ao=25°  and  ai=5°.  Now,  the  perturbations  of  the  sweep  angle  with  the  same 
frequency  as  pitching,  but  with  a  phase  angle  are  considered,  i.e. 
a=ao+ai  sin  ©t 
A=Ao+Ai  sin  (©t+O) 

The  variation  of  breakdown  location  for  three  different  values  of  the  phase  angle  O  is 
shown  in  Figure  30.  The  phase  angle  <I)  was  varied  over  a  range  of  360°,  and  the  maximum 
and  minimum  values  of  breakdown  location  during  a  cycle  was  documented.  The  results 
are  shown  in  Figure  31.  It  appears  that  for  0=0°,  the  breakdown  location  is  approximately 
constant  and  the  average  location  has  moved  downstream.  (The  dashed  lines  indicate  the 
maximum  and  minimum  location  for  the  pitching  motion  only  at  the  average  sweep  angle 
Aq).  The  perturbations  of  sweep  angle  not  only  decrease  the  amplitude  of  the  variations  of 
breakdown  location,  but  also  move  the  average  location  downstream. 

The  results  for  a  different  set  of  parameters  (Ao=66°,  Aj=3°,  ao=20°  and  ai=5°)  are 
shown  in  Figure  32.  In  this  case,  the  amplitude  of  the  variations  decreases  considerably  at 
an  optimum  phase  angle  around  35°.  However,  the  average  location  does  not  change  much. 
Experiments  for  a  different  range  of  angle  of  attack  and  sweep  angle,  as  well  as  reduced 
frequency,  showed  that  the  optimum  phase  angle  is  somewhere  between  0°  and  35°.  A 
simple  explanation  of  the  optimum  phase  angle,  based  on  the  time-lag  response  of 
breakdown  location  (with  respect  to  the  variations  in  a(t)  and  A(t)),  is  possible.  As  a  first 
approximation,  the  response  of  vortex  breakdown  location  to  harmonic  variations  of  angle 
of  attack  (a=ao+ai  sin  ©t,  for  constant  sweep  angle  A=Ao)  will  be  in  the  form  of  Xbd=Xo+Xi 


sin  (cot  +7i-(j)J,  where  <j)„  is  the  phase  lag  with  respect  to  the  quasi-steady  case.  Since,  in 
the  quasi-steady  case,  the  vortex  breakdown  location  moves  toward  apex  with  increasing 
angle  of  attack,  there  is  a  phase  angle  of  n  in  the  expression  of  x^j.  Similarly,  the  response 
of  vortex  breakdown  location  to  harmonic  variations  of  sweep  angle  (A=Ao+Ai  sin  (cot+O), 
for  constant  angle  of  attack  a=ao)  will  be  in  the  form  of  Xb(j=Xo+Xi  sin  (0t+O-<|)^),  where 
<j)y^  is  the  phase  lag  with  respect  to  the  quasi-steady  case.  Since  it  is  desired  to  cancel  the 
variations  of  breakdown  location  due  to  a(t)  with  the  variations  due  to  A(t),  the  arguments 
of  these  two  responses  should  be  out  of  phase,  i.e.  ((0t+O-(l)4)+rt=((0t  +7t-(|)„),  which  yields 
<I)=(|)^-(|)a.  For  the  case  shown  in  Figure  32,  the  phase  lags  (j)/^  and  (|)oi  were  found  by 
plotting  the  time  history  of  breakdown  location,  which  yielded  (|)a==85°  and  <|)a=40°. 
Therefore,  O  is  predicted  as  0=(t)A-(l)(j=45°,  which  is  close  to  the  experimental  results 
shown  in  Figure  32.  Although  this  simple  model  provided  a  reasonable  prediction  for  the 
optimum  phase  angle  between  a(t)  and  A(t),  detailed  flow  field  information  is  necessary  for 
further  understanding. 

In  summary,  harmonic  variations  of  sweep  angle  were  studied  for  a  delta  wing  with 
variable  sweep.  The  variation  of  vortex  breakdown  location  showed  hysteresis  loops  and 
phase  lag  which  depend  on  the  reduced  frequency.  However,  measurements  of  the  phase- 
averaged  axial  velocity  in  a  cross  flow  plane  upstream  of  breakdown  location  showed  no 
major  vmsteady  effects  in  the  development  of  the  flow  field.  The  larger  time  lag  of  vortex 
breakdown  location  cannot  be  explained  with  the  variations  in  the  upstream  conditions 
since  the  location  and  strength  of  the  vortex  exhibit  smaller  phase  lags.  This  suggests  that 
the  external  pressure  gradient  plays  a  major  role  in  the  dynamic  response  of  breakdown 
location.  Use  of  variable  sweep  for  breakdown  control  was  demonstrated  for  a  pitching 
delta  wing.  Oscillations  of  sweep  angle  with  the  same  frequency  as  pitching,  but  with  a 
phase  angle  were  considered.  It  was  shown  that,  at  an  optimum  phase  angle,  the 
perturbations  of  sweep  angle  decrease  the  amplitude  of  the  variations  of  breakdown 
location.  In  some  cases,  it  was  also  observed  that  the  average  breakdown  location  moved 
downstream. 


CONCLUSIONS 


Coherent  pressure  fluctuations  observed  on  delta  wings  are  due  to  the  first  helical 
mode  instability  of  the  swirling  flows.  Characteristics  of  this  instability  was  studied  in 
detail. 

The  fluctuations  of  vortex  breakdown  location  occur  at  much  lower  frequencies  than 
the  frequency  of  the  hydrodynamic  instability  of  the  flow  in  the  wake  of  breakdown.  It  is 
thus  suggested  that  the  helical  mode  instability  does  not  influence  the  unsteady  nature  of 
breakdown  location.  The  source  of  the  fluctuations  was  identified  as  the  interaction 
between  the  opposite  leading-edge  vortices. 

A  feedback  control  loop  was  used  to  control  vortex  breakdown  location.  For  this 
purpose,  the  amplitude  of  pressure  fluctuations  was  chosen  as  the  control  variable.  It  was 
demonstrated  that  it  is  feasible  to  use  a  variable  sweep  angle  mechanism  as  a  means  of 
controlling  the  vortex  breakdown  location  by  influencing  the  circulation  of  the  leading  edge 
vortex. 

Unsteady  pressure  measurements  were  carried  out  on  a  pitching  delta  wing  in  order 
to  study  the  effects  of  pressure  gradient  on  the  phase  lag  between  the  wing  motion  and  the 
movement  of  breakdown  location.  The  results  indicate  that  the  observed  time  lag  of 
breakdown  location  is  strongly  linked  to  the  external  pressure  gradient  generated  by  the 
wing. 

Effect  of  leading-edge  devices  on  vortex  breakdown  was  investigated.  Stationary  as 
well  as  oscillating  leading-edge  flaps  were  shown  to  alter  the  structure  of  the  leading-edge 
vortices,  and  consequently  the  location  of  vortex  breakdown.  Also,  periodic  variations  of 
the  sweep  angle  were  used  to  control  the  breakdown  location  for  a  pitching  delta  wing. 
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Figure  4.  Experimental  setup 
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Figure  5.  Location  of  vortex  breakdown  and  pressure  fluctuations  near  trailing-edge 


(x/c=0.93)  for  different  angles  of  attack.  Length  of  time  record  is  approximately  lOc/Uoc. 


Phase  lag 


ISO- 


c 


*-*.*»«  Qlpha=3l.3  deg 
alpha=33.3  deg 
‘iSPqp  alpha=35.3  deg 


0.20  0.40  0.60  0.80 

y/s 


Figure  6.  Variation  of  phase  angle  between  pressure  siganals  at  two  locations  in  the 
span  wise  direction. 


Figure  9.  Cross-spectrum  of  flucmation  of  left  and  right  breakdowns,  A=70‘’  and  a=3T. 
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Figure  10.  Spectrum  of  unsteady  flow  phenomena  over  delta  wings. 
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Figure  14.  Pressure  fluctuation  and  time  history  of  angle  of  attack  for  a  periodically 
pitching  motion.  Length  of  time  record  is  838ms;  x/c=0.55,  y/s=0.5,  k=0.35. 


Figure  15.  Variation  of  standard  deviation  as  a  function  of  angle  of  attack 
at  different  streamwise  locations. 
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Figure  16.  Mean  pressure  coefficient  in  spanwise  direction. 
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Figure  17.  Variation  of  phase-averaged  pressure  along  centerline. 
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Figure  21.  Constant  contours  of  normalized  axial  velocity  for  5=0"  (top,  a=20"),  5=60' 
(middle,  a=25"),  5=100°  (bottom,  a=20°). 
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Figure  22:  Variation  of  maximum  swirl  velocity  as  a  function  of  flap  angle  for  a=20°  and 
a=25°.  The  arrows  indicate  the  flap  angle  at  which  breakdown  is  at  the  trailing-edge  of 
the  wing. 


Figure  23:  Variation  of  vortex  breakdown  location  for  harmonic  variations  of  flap  angle 
(k=0.4),  and  for  static  flaps  (a=30°). 
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Figure  27:  Contours  of  st 


on  and  rms  velocity  (in  the  static  case). 


Figure  29:  Variation  of  breakdown  location  for  combined  motion  of  pitching  and  variable 
sweep  for  different  values  of  phase  angle. 
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Figure  30.  Vatiation  of  breakdown  location  for  combined  motion  of  pitching  and  variable 

sweep  for  different  values  of  phase  angle. 
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Figure  31:  Maximum  and  minimum  locations  of  breakdown  as  a  function  of  phase  angle. 
The  dashed  line  is  for  pitching  motion  only  for  the  average  sweep  angle. 


Figure  32:  Maximum  and  minimum  locations  of  breakdown  as  a  function  of  phase  angle. 
The  dashed  line  is  for  pitching  motion  only  for  the  average  sweep  angle. 
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Unsteady  Flow  Phenomena  over  Delta  Wings  at 
High  Angle  of  Attack 


Ismet  Gursul* 

University  of  Cincinnati,  Cincinnati,  Ohio  45221 


Experiments  show  that  coherent  pressure  fluctuations  observed  on  delta  wings  are  due  to  the  helical  mo  e 
instability  of  the  vortex  breakdown  Howfield.  No  dominant  frequency  in  the  spectra  of  pressure  Huctuations  on 
the  wing  surface  was  observed  after  the  breakdown  reached  the  apex  of  the  wing,  aithough  the  vortex  shedding 
could  be  detected  in  the  wake.  Measurements  of  pressure  fluctuations  at  different  streamwise  locations  on  the 
wing  suggest  that  the  dimensionless  frequency /x/f/^  is  nearly  constant  for  a  given  pometry,  which  implies 
increasing  wavelength  in  the  streamwise  direction.  For  different  wings,  this  nondimensional  frequency  is  shown 
to  be  a  function  of  nondimensional  circulation  r/i;„  x  only.  Both  the  wavelength  of  the  disturbances  and  the  ^re 
radius  increase  with  the  nondimensional  circulation  at  a  fixed  streamwise  location.  The  wavelength  norma  ize 
by  the  core  radius  is  around  3-4,  which  is  much  smaller  than  the  predictions  for  the  Q  vortex. 


Nomenclature 

c  =  root  chord 
/  =  frequency 

k  =  axial  wave  number 

n  =  azimuthal  wave  number 

p  =  pressure 

q  =  parameter  for  Q  vortex 
Re  =  Reynolds  number 
r  =  radial  distance  from  vortex  axis 
Tq  =  characteristic  core  size 
s  =  local  semispan 
=  freestream  velocity 
fj^  =  phase  speed 
U,  =  axial  velocity  excess  or  deficit 
V  =  azimuthal  velocity 
W  =  axial  velocity 

X  =chordwise  distance  from  wing  apex 
y  =  spanwise  distance  from  wing  root 
2  =  distance  above  wing  surface 

a  =  angle  of  attack 

r  =  leading-edge  vortex  circulation 
0  =  phase  angle 

A  =  sweep  angle 

X  =  wavelength 

p  =  density 

<I>  =  phase  angle 

^  azimuthal  angle 
CD  =  radial  frequency 

Introduction 

LTHOUGH  interest  in  high  angle-of-attack  aerodynamics  has 
increased,  little  is  known  about  unsteady  flow  in  vortex 
breakdown/postbreakdown  flowfields.  Sharply  increased  fluctua¬ 
tions  in  the  normal-force  coefficient  for  delta  wings  were  observed 
when  vortex  breakdown  moved  over  the  wing.^  Measurements  of 
surface  pressure  fluctuations  in  the  wake  of  breakdown  showed 
coherent  oscillations.^’^  Periodic  oscillations  were  also  observed  in 
other  swirling  flows  after  breakdown  occurred.^^  Garg  and  Leibo- 
vich  carried  out  LDV  measurements  in  the  wakes  of  breakdown  in 
tubes  and  observed  coherent  oscillations.  They  suggested  that  the 
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measured  frequencies  correspond  to  the  theoretical  predictions 
found  in  Ref.  8  for  the  first  helical  mode  (n  =  - 1)  of  the  time-aver¬ 
aged  mean  flow  profiles,  assuming  that  the  oscillations  are  the  dis¬ 
turbances  with  the  maximum  growth  rate,  [The  disturbances  are 
represented  as  exp{i(Ax  +  -cor)),  where  o)  is  the  frequency,  k 

the  wave  number  in  the  axial  direction,  and  n  the  wave  number  (an 
integer)  in  the  angular  direction].  Thus,  one  possible  source  of  the 
observed  oscillations  over  delta  wings  is  the  hydrodynamic  insta¬ 
bility  of  breakdown  wake. 

On  the  other  hand,  based  on  the  velocity  measurements  in  the 
wake  of  a  delta  wing,  vortex  shedding  was  shown  to  occur.’  It  Was 
suggested  that  the  alternate  shedding  at  large  ^gle  of  at^k 
should  induce  asymmetry  on  the  pressure  distribution  on  the  wing. 
The  smallest  angle  of  attack  for  which  the  periodic  motion  was 
detected  was  around  35  deg  for  the  particular  delta  wing  used. 
Although  the  authors  do  not  report  the  location  of  vortex  bre^- 
down,  it  is  estimated  to  be  over  the  wing  and  close  to  the  trailing 
edge.*®’"  This  suggests  that  vortex  shedding  may  start  when  the 
breakdown  moves  over  the  wing.  Although  the  time-averaged 
velocity  measurements  (for  example,  see  Hummel*^)  reveal  that 
the  swirling  flow  persists  in  the  wake  of  breakdown,  little  is  known 
about  the  instantaneous  flow  structure.  Since  the  axial  convection 
of  vorticity  along  the  core  is  reduced  after  the  breakdown,  vortex 
shedding  might  start  from  the  wing. 

The  purpose  of  this  study  is  to  investigate  the  characteristics  of 
unsteady  flow  over  delta  wings  of  different  sweep  angle.  The 
nature  of  coherent  oscillations,  their  origin  (helical  mode  instabil¬ 
ity  vs  vortex  shedding),  and  their  influence  on  unsteady  loading 
will  be  discussed  with  the  aid  of  measurements  of  pressure  and 
velocity  as  well  as  flow  visualization. 

Experimental  Setup 

The  experiments  were  carried  out  in  an  open-circuit  wind  tunnel 
having  a  cross  section  of  305  mm  X  305  mm.  The  flow  entered 
through  a  9.5:1  contraction  section  with  honeycomb  and  screens, 
passed  through  the  test  section,  and  was  exhausted  by  an  axial  fan. 
The  turbulence  intensity  in  the  test  section  was  about  0.25%.  The 
freestream  velocity  spectrum  was  free  from  any  sh^  peaks.  The 
angle  of  attack  of  the  wings  could  be  varied  continuously  from 
outside  the  tunnel. 

Pressure  and  velocity  measurements  were  made  for  four  delta 
wings  with  sweep  angles  A  =  60,  65,  70,  and  75  deg.  The  chord 
lengths  ranged  from  120  mm  to  150  mm,  and  the  models  were  con¬ 
structed  of  9.5-mm-thick  PVC.  The  lee  surfaces  were  flat  whereas 
the  leading  edges  were  beveled  at  45  deg  on  the  windward  side. 
The  models  were  sting  mounted.  The  measurements  were  done  at 
Reynolds  numbers  based  on  the  chord  length  equal  to  25,000- 
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Fig.  1  Spectrum  of  pressure  fluctuations  at  x/c  =  0.89,  y/s  =  0.42,  a  = 
353  deg,  A  =  70  deg;  vertical  scale  is  logarithmic. 


Fig.  2  Phase  delay  between  the  pressure  fluctuations  measured  at 
points  1  and  2;  point  i  is  fixed,  while  the  location  of  point  2  is  varied 
(top).  Variation  of  the  amplitude  of  pressure  fluctuations  (bottom),  A  = 
70  deg. 

lOO.CXX),  although  most  of  the  results  reported  here  correspond  to 
Re  =  100,(X)0.  At  the  largest  angle  of  attack,  the  blockage  ratio 
was  0.094.  A  piezoelectric  pressure  transducer  (PCB,  model 
103A)  was  used  to  measure  the  fluctuating  pressure  on  the  wing 
surface.  The  transducer  has  been  flush  mounted  at  the  surface.  The 
transducer  and  its  leads  were  housed  in  the  grooves  running  either 
along  a  ray  from  the  apex  to  the  trailing  edge  (at  40-60%  of  the 
local  semispan,  depending  on  the  sweep  angle)  or  along  the  span- 
wise  direction.  The  pressure  sensing  part  of  the  transducers  had  a 
diameter  of  2.54  mm.  Two  pressure  transducers  were  used  for 


simultaneous  measurements  in  the  spanwise  and  streamwise  direc¬ 
tions.  A  hot-wire  probe  was  used  to  obtain  the  spectra  of  the  veloc¬ 
ity  fluctuations  over  the  wing  and  in  the  wake.  Pressure  and  veloc¬ 
ity  signals  were  processed  by  a  two-channel  signal  analyzer  (HP 
35660A).  Smoke-wire  visualization  by  a  0.1-mm-diam  stainless 
steel  wire,  and  the  smoke  injected  near  the  apex  provided  informa¬ 
tion  on  the  flow  structure  and  vortex  breakdown  location.  A  light 
source  consisting  of  two  600-W  lamps  was  used  to  illuminate  the 
flowfield.  In  addition  to  the  still  pictures  taken,  flow  visualization 
was  videotaped  at  low  speeds.  The  measurement  uncertainty  for 
dimensionless  pressure  coefficient  ranged  from  10%  to  0.6%, 
whereas  the  uncertainty  for  frequency  and  breakdown  location 
were  estimated  as  1.4%  and  3%,  respectively. 

Results 

When  the  vortex  breakdown  moved  over  the  wing,  the  pressure 
fluctuations  exhibited  coherent  oscillations.  An  example  of  pres¬ 
sure  fluctuation  spectrum  for  sweep  angle  A  =  70  deg,  at  an  angle 
of  attack  a  =  35.3  deg  is  shown  in  Fig.  1.  The  measurement  was 
taken  at  x/c  =  0.89,  y/s  =  0.42  (where  s  is  the  local  semispan)  and 
the  vortex  breakdown  location  was  -  0.3.  A  sharp  peak  in  the 
spectra  was  always  observed  in  the  wake  of  breakdown  location, 
regardless  of  angle  of  attack  and  sweep  angle  for  all  of  the  delta 
wings  tested.  In  the  Reynolds  number  range  investigated,  the 
dimensionless  dominant  frequency/c/f/«  seemed  to  be  independent 


Fig.  3  Smoke  visualization  showing  helical  instability  in  the  break¬ 
down  region  (a  =  35  deg,  A  =  70  deg),  Re  =  15,000. 
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of  the  Reynolds  number.  To  understand  the  nature  of  these  oscilla¬ 
tions,  measurements  were  taken  with  two  pressure  transducers 
located  in  the  spanwise  direction.  Cross-spectral  analysis  provided 
the  phase  angle  between  the  pressure  fluctuations  measured  with 
the  two  transducers.  If  the  pressure  fluctuations  are  represented  as 
cos  (CD/)  and  cos  (cor  —  9)  at  points  I  and  2,  respectively,  the  phase 
angle  0  is  shown  in  Fig.  2  together  with  the  amplitude  (Cp  = 
distribution  in  the  spanwise  direction.  In  measuring 
the  phase  angle,  the  location  of  point  1  was  fixed  while  the  legation 
of  point  2  was  varied.  Increase  in  phase  delay  toward  the  leading 
edge  is  consistent  with  the  direction  of  rotation  in  the  vortex.  'Hie 
vortex  shedding  from  the  leading  edge  cannot  cause  this  v^ation 
of  the  phase  delay.  Instead,  the  trend  would  likely  be  opposite  (i.e., 
decreasing  phase  lag  toward  the  leading  edge).  Thus,  the  evidence 
suggests  that  the  oscillations  may  be  due  to  the  helical  mode  insta¬ 
bility.  Indeed,  the  flow  visualization  obtained  by  releasing  smoke 
near  the  apex  shows  this  type  of  instability  (Fig.  3).  In  addition,  the 
azimuthal  wave  number  was  estimated  from  the  phase  angle 
between  the  signals  from  two  hot-wire  probes  located  at  y/s  =  0.5 
and  z/s  =  0.25  and  y/s  =  0.5  and  z/s  =  0.75.  With  respect  to  the  vor¬ 
tex  axis,  these  locations  are  approximately  180  deg  apart.  The 
probes  were  oriented  so  as  to  be  most  sensitive  to  tangential  and 
axial  velocities.  The  phase  angle  is  shown  in  Fig.  4  for  several  val¬ 
ues  of  angle  of  attack,  which  suggests  that  the  azimuthal  wave 
number  is /I  =  1. 

These  results  are  consistent  with  the  predictions  of  the  linear 
stability  theory.®  Lessen  et  al.  carried  out  their  calculations  for  a 
jet-like  axial  velocity  profile  and  showed  that  the  flow  is  most  sen¬ 
sitive  to  the  disturbances  with  negative  azimuthal  wave  numbers 
(n  <  0).  Translation  and  inversion  of  the  axial  velocity  profile  does 
not  affect  the  growth  rate  within  the  temporal  analysis.  However, 
the  azimuthal  wave  number  changes  its  sign.^^  Therefore,  for  a 
wake-like  profile  such  as  found  in  breakdown  region,  the  most 
unstable  modes  are  the  ones  with  positive  azimuthal  wave  num¬ 
bers  {n  >  0),  which  represent  disturbances  rotating  in  the  same 
direction  as  the  vortex.  This  is  in  agreement  with  the  findings  in 
Fig.  2.  The  theory  cannot  predict  which  azimuthal  wave  number 
will  be  amplified.  However,  experiments  in  a  tube'^  and  tip  vortex 
showed  that  n=\.  This  is  also  in  agreement  with  the  present 
results  (see  Fig.  4).  The  Ini  =  I  mode  is  necessary  for  a  streakline 
released  on  the  axis  to  take  a  helical  shape,'^  since  the  radial  veloc¬ 
ity  component  must  be  zero  on  the  axis  for  all  modes  except  Ini  =  1 . 
In  the  breakdown  region  shown  in  Fig.  3,  the  large-scale  helical 
shape  of  smoke,  which  is  originated  from  the  axis,  is  an  evidence 
of  the  helical  mode  instability.  Since  the  disturbances  with  n  =  1 
will  have  a  phase  function  (fcx  +  ()>  -  cor),  one  expects  velocity/vor- 
ticity  to  be  periodical  in  the  x  direction  and  be  antisymmetric  with 
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——  Dominant  frequency 

»  «  *  »  *  Vortex  breakdown  location 


65  deg. 


Fig.  6  Vortex  shedding  at  large  angle  of  attack  (a  =  60  deg),  A  =  60 
deg  and  A  =  75  deg.  Re  =  15,000, 


respect  to  the  vortex  axis.  Indeed,  the  instantaneous  azimuthal  vor- 
ticity  in  a  plane  through  the  vortex  axis'*  shows  that  vorticity  con¬ 
centrations  form  the  two  rows  of  periodic  structures  located  on 
either  side  of  the  axis  with  an  antisymmetric  configuration,  resem¬ 
bling  the  well-known  Karman  vortex  street. 

With  increasing  angle  of  attack,  the  pressure  spectra  became 
free  from  any  sharp  peaks.  The  domain  over  which  a  sharp  peak 
was  observed  in  the  spectrum  for  the  wing  with  sweep  angle  A  = 
65  deg  is  shown  in  Fig.  5.  A  dominant  frequency  was  obrerved  at 
the  measurement  stations  downstream  of  breakdown  location,  until 
the  breakdown  reached  the  apex  of  the  wing.  When  the  angle  of 
attack  is  larger  than  the  one  that  corresponds  to  the  breakdown 
location  at  apex,  no  oscillations  were  observed  in  pressure  fluctua¬ 
tions.  The  measurements  for  other  wings  gave  results  similar  to 
Fig.  5.  Thus,  it  is  concluded  that,  until  the  breakdown  reaches  the 
apex  of  the  wing,  swirling  flow  persists  over  the  wing  and  hydro- 
dynamic  instability  of  this  flow  is  the  only  source  of  the  unsteady 
loading  on  the  wing.  For  larger  angles  of  attack,  the  shear  layer 
separated  from  the  leading  edge  will  not  be  able  to  form  a  swirling 
flow  with  axial  motion.  Thus,  at  large  angle  of  attack,  the  vortex 
shedding  from  the  wing  starts  as  shown  in  Fig.  6.  Whereas  the 
shedding  seems  symmetric  for  the  large  aspect  ratio  wing,  alter¬ 
nate  shedding  becomes  dominant  for  the  low  aspect  ratio  wing.  A 
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Fig.  7  Variation  of  dimensionless  dominant  frequency  on  the  wing 
{xfc  =  0.89)  and  at  one  chord  length  downstream  of  trailing  edge  as  a 
function  of  angle  of  attack,  A  =  60  deg. 


x/ c 


Fig.  8  Variation  of  dimensionless  frequency  fciV^  as  a  function  of 
streamwise  distance. 


Measurements  of  pressure  at  different  streamwise  locations  on 
the  wing  upper  surface  showed  that  the  dominant  frequency  is  not 
constant,  but  depends  on  the  distance  from  the  apex.  Examples  of 
variation  of  dimensionless  frequency  fclU„  as  a  function  of  the 
streamwise  distance  for  different  wings  is  shown  in  Fig.  8.  The 
scatter  in  the  data  is  believed  to  be  due  to  fluctuations  in  the  break¬ 
down  location.  It  is  well  known  that  vortex  breakdown  location  ex¬ 
hibits  irregular  variations.  The  dominant  frequency  decreases  with 
increasing  distance.  This  was  also  noticed  by  Roos  and  Kegelman^ 
(for  a  wing  with  A  =  70  deg  and  a  =  33  deg);  moreover,  it  was 
shown  that  if  the  frequency  is  normalized  with  the  distance  jc,  the 
dimensionless  frequency  ficlU^  is  nearly  constant.  This  was  attrib¬ 
uted  to  disturbances  rotating  in  the  vortex  core  whose  radius  grows 
linearly  with  jc.  The  variation  of  fxlU^  corresponding  to  the  data  in 
Fig.  8  is  shown  in  Fig.  9.  It  seems  that  fxlV^  is  nearly  constant  for 
a  given  A  and  a  for  all  delta  wings  tested.  Since  the  disturbances 
are  shown  to  be  helical  waves,  this  implies  that  the  wavelength  in¬ 
creases  with  X,  Indeed,  the  instantaneous  azimuthal  vorticity  in  a 
plane  through  the  vortex  axis,  which  was  recently  obtained  by  a 
particle  image  velocimetiy  technique,'^  clearly  shows  vorticity 
concentrations  with  increasing  wavelength  between  them.  The  in¬ 
stantaneous  azimuthal  vorticity  concentrations  seem  to  be  located 
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Fig.  9  Variation  of  dimensionless  frequency  fxiV^  as  a  function  of 
streamwise  distance. 
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signal  from  a  hot-wire  probe  located  at  one  chord  length  from  the 
trailing  edge  was  used  to  monitor  unsteady  flow  in  the  wake.  Dom¬ 
inant  frequency  detected  in  the  wake  was  compared  with  the  dom¬ 
inant  frequency  detected  by  the  pressure  transducer  located  near 
the  trailing  edge  {xlc  =  0.89)  in  Fig.  7.  Until  the  breakdown  loca¬ 
tion  reaches  the  apex  (a=  30  deg  for  the  wing  with  sweep  angle  A 
=  60  deg),  both  signals  have  the  dominant  frequency  of  the  helical 
mode  instability,  although  the  frequencies  differ  slightly.  This  may 
be  due  to  the  streamwise  gradients  of  the  mean  flow  in  the  wake. 
Since  no  coherent  pressure  oscillations  on  the  wing  due  to  vortex 
shedding  was  detected  (see  also  Fig.  5),  this  suggests  that  unsteady 
loading  on  the  wing  due  to  vortex  shedding  is  negligible.  After  the 
helical  mode  instability  disappears,  vortex  shedding  frequency  is 
detected  in  the  wake.  Although  the  relationship  between  the  domi¬ 
nant  frequency  in  the  wake  and  angle  of  attack  seems  a  continuous 
curve,  it  covers  two  different  flow  regimes.  Indeed,  in  the  vortex 
shedding  regime,  the  frequency  is  pretty  much  constant.  Similar 
trend  was  observed  for  a  wing  with  A  =  76  deg  in  other  experi¬ 
ments.® 
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Fig.  10  Variation  of  dimensionless  frequency  fxIU^  as  a  function  of 
angle  of  attack  for  different  values  of  sweep  angle. 
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angle  of  attack  and  sweep  angle  was  estimated  by  using  a  slender 
wing  theory.*^  Figure  11  shows  that  all  of  the  curves  cojlapse 
except  for  A  =  60  deg.  Since  the  method  used  to  cdculate  circula¬ 
tion  is  a  good  approximation  for  low  aspect  ratio  wings,  the  dis¬ 
crepancy  for  A  =  60  deg  is  believed  to  be  due  to  the  incapability  of 
the  method  to  estimate  the  circulation. 

The  wavelength  of  the  disturbance  was  measured  by  two  pres¬ 
sure  transducers  located  in  the  streamwise  direction.  From  the 
phase  measurements  obtained  by  the  cross-spectral  analysis,  the 
wavelength  was  calculated  as  X  =  27i/(3<l>/3x),  where  <1>  is  mea¬ 
sured  in  radians.  Typical  distance  between  the  transducers  was  Ax/ 
c  =  0.066,  corresponding  to  Ax/X=  0.13.  The  midpoint  between  the 
transducers  was  located  near  the  trailing  edge  {x/c  =  0.85).  Based 
on  the  previous  results,  the  wavelength  is  expected  to  be  linearly 
changing  with  x  in  the  conical  mean  flowfield.  The  variation  of  X/x 
as  a  function  of  angle  of  attack  and  dimensionless  circulation  is 
shown  in  Fig.  12.  In  addition,  the  variation  of  the  phase  speed  Uc  = 
fX  is  shown  in  Fig.  13.  The  quantity  X/x  seems  to  be  a  linear  func¬ 
tion  of  the  dimensionless  circulation.  A  linear  curve  fit  provided 
the  following  relation: 


(2) 


Since  the  normalized  convection  speed  seems  to  be  inde¬ 

pendent  of  the  dimensionless  circulation,  it  can  be  shown  that 

^  ^  ^  ^ 

X  fx  fx  ifx/uj 

and  using  £q.  (2) 

X  1  r 
^  ^  ifx/UJ  “  U^x 

Hence, 


Indeed,  a  curve  fit  to  the  data  shown  in  Fig.  11  gives  the  following 
relationship: 


^  =  1.038 
£/„ 


Fig.  15  Variation  of  X/r^  as  a  function  of  parameter  q. 


Lessen  et  al.®  calculated  the  most  unstable  wavelength  for  the 
first  helical  mode  of  the  following  time-averaged  mean  flow: 

V/U^  =  -4-n-exp(-rVro)] 


2  2 

W/U^  =  exp(-r  /g 

where  the  circulation  is  T  =  InqrQU,.  This  is  known  as  Q  vortex 
and  the  parameter  q  is  directly  proportional  to  the  circulation  and 
inversely  proportional  to  the  axial  velocity  excess  or  deficit 
The  maximum  swirl  speed  occurs  at  r  =  1.122  Tq.  Therefore,  the 
parameter  Tq  can  be  taken  as  the  characteristic  core  size.  The  most 
unstable  wavelength  X/rQ  for  the  first  helical  mode  was  calculated 
by  Lessen  et  al.®  for  several  values  of  the  parameter  q.  The  core  ra¬ 
dius  To  is  expected  to  be  a  linear  function  of  x  in  the  conical  flow 
over  delta  wings.  From  a  limited  number  of  mean  velocity  mea¬ 
surements,^'*®-’^  To  lx  was  estimated  and  is  shown  in  Fig.  14.  The 
data  suggests  that  the  core  radius  increases  with  the  increasing 
swirl  level,  which  was  predicted  for  a  trailing  vortex^*^  and  swirling 
flows  in  tubes,^*  Although  the  data  is  very  limited,  there  is  a  trend 
of  a  linear  relationship  between  r^lx  and  T lU^x,  This  together  with 
Eq.  (2)  suggests  that  the  normalized  wavelength  X/kq  is  roughly 
constant  as  is  expected  in  a  conical  flowfield.  The  variation  of  Tl/fq 
with  the  estimated  parameter  q  is  shown  in  Fig.  15,  which  indi¬ 
cates  that  the  disturbances  have  a  much  smaller  wavelength  than 
the  predictions  for  the  Q  vortex.  This  may  be  due  to  slowly  varying 
mean  flow  in  the  streamwise  direction.  Garg  and  Leibovich** 
showed  that  the  measured  frequencies  in  vortex-tube  experiments 
agree  very  well  with  the  temporal  stability  predictions  for  the  Q 
vortex,  assuming  that  the  disturbances  have  the  wavelength  and 
phase  speed  corresponding  to  the  maximum  growth  rate.  On  the 
other  hand,  the  experimental  result  for  a  trailing  vortex’^  (for  which 
the  streamwise  gradient  is  also  negligible)  is  closer  to  the  results 
for  leading-edge  vortex  over  delta  wings.  The  discrepancy  be¬ 
tween  the  predictions  for  Q  vortex  and  the  present  experimental  re¬ 
sults  remains  to  be  explained. 

Conclusions 

Experiments  showed  that  coherent  pressure  fluctuations  are 
observed  on  delta  wings  as  long  as  vortex  breakdown  is  over  the 
wing.  Cross-spectral  analysis  of  two-point  pressure  and  velocity 
measurements  as  well  as  flow  visualization  confirm  that  these 
oscillations  are  in  the  form  of  helical  waves  with  azimuthal  wave 
number  n  =  1,  which  is  the  result  of  the  instability  of  swirling 
breakdown  wake  flow.  The  source  of  unsteady  loading  on  the  wing 
is  due  to  this  helical  mode  instability  as  opposed  to  vortex  shed¬ 
ding  which  starts  after  the  breakdown  reaches  the  apex  of  the 
wing.  Vortex  shedding  takes  the  form  of  symmetric  or  alternate 
depending  on  the  aspect  ratio  of  the  wing.  The  influence  of  vortex 
shedding  on  unsteady  pressure  fluctuations  is  negligible  and  could 
not  be  detected  in  the  experiments,  although  the  velocity  fluctua¬ 
tions  in  the  wake  show  coherent  oscillations. 

Measurements  of  unsteady  pressure  at  different  streamwise 
locations  on  the  wing  surface  showed  that  the  dimensionless  fre¬ 
quency  fic/U^  is  nearly  constant  for  a  fixed  angle  of  attack  and 
sweep  angle.  Based  on  the  results  from  this  work  and  other  experi¬ 
ments,  it  is  suggested  that  the  wavelength  of  the  disturbances 
increases  linearly  with  jc  in  the  conical  mean  flowfield,  resulting  in 
constant  fx/U„.  The  nondimensional  frequency  is  found  to  be  a 
function  of  nondimensional  circulation  T/U^x  only,  where  the  cir¬ 
culation  was  calculated  by  assuming  conical  flow  and  using  a  slen- 
der-wing  theory. 

Measurements  of  the  wavelength  and  phase  speed  indicates  that 
X/x  increases  with  increasing  TIU^x,  whereas  the  phase  speed  is 
roughly  constant  -  0.6.  The  core  radius  Fq  which  is  approx¬ 
imately  equal  to  the  radius  at  which  the  swirl  speed  becomes  max¬ 
imum,  also  increases  with  increasing  T/U^x.  The  normalized 
wavelength  XJfq  seems  to  be  around  3-4,  which  is  much  smaller 
than  the  predictions  for  the  Q  vortex. 
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Finally,  it  should  be  noted  that  the  existence  of  helical  mode  in¬ 
stability  does  not  have  any  implications  with  regard  to  the  break¬ 
down  mechanism.  The  flows  upstream  of  breakdown  location  are 
known  to  be  stable  to  disturbances.  However,  whether  the  instabil¬ 
ities  in  downstream  region  promote  breakdown  is  not  known. 

Acknowledgments 

Part  of  this  work  was  supported  by  the  Air  Force  Office  of  Sci¬ 
entific  Research.  The  author  would  like  to  thank  F.  Roos  for  mak¬ 
ing  some  of  his  data  available. 

References 

^Eamshaw,  P.  B.,  and  Lawford,  J.  A.,  “Low  Speed  Wind  Tunnel  Exper¬ 
iments  on  a  Series  of  Sharp-edged  Delta  Wings,”  Aeronautical  Research 
Council,  R  &  M  3424,  Aug.  1964. 

^Mabey,  D,  B.,  “Beyond  the  Buffet  Boundary.”  Aeronautical  Journal, 
April  1973,  pp.  201-215. 

^Roos,  F.  W.,  and  Kegelman,  J.  T.,  “Recent  Explorations  of  Leading- 
edge  Vortex  Rowfields,”  NASA  High-Angle-of-Attack  Technology  Con¬ 
ference,  NASA  Langley  Research  Center,  Hampton,  VA,  Oct.  30-Nov.  1, 
1990. 

'^Garg,  A.  K.,  and  Leibovich,  S.,  “Spectral  Characteristics  of  Vortex 
Breakdown  Flowfields,”  Physics  of  Fluids,  Vol.  22,  No.  11,  1979,  pp. 
2053-2064. 

^Chanaud,  R.  C.,  “Observations  of  Oscillatory  Motion  in  Certain  Swirl¬ 
ing  Flows,”  ioMrna/  of  Fluid  Mechanics,  Vol.  21,  1965,  pp.  111-127. 

^Cassidy,  J.  J.,  and  Falvey,  H.  T.,  “Observations  of  Unsteady  Flow  Aris¬ 
ing  after  Vortex  Breakdown,”  Journal  of  Fluid  Mechanics,  Vol.  41,  1970, 
pp.  727-736. 

‘^Singh,  P.  I.,  and  Uberoi,  M.  S.,  “Experiments  on  Vortex  Stability,” 
Physics  of  Fluids,  Vol.  19,  No.  12, 1976,  pp.  1858-1863. 

^Lessen,  M.,  Singh,  P.  J.,  and  Paillet,  F.,  “The  Stability  of  a  Trailing 
Line  Vortex,  Part  1,  Inviscid  Theory,”  Journal  of  Fluid  Mechanics,  Vol. 


63, 1974,  pp,  753-763. 

^ediniotis,  O.  K.,  Stapountzis,  H.,  and  Telionis,  D.  P.,  “Vortex  Shed¬ 
ding  over  Delta  Wings,”  AJAA  Journal,  Vol.  28,  No.  5, 1990,  pp.  944-946. 

‘^entz,  W.  H.,  and  Kohlman,  D.  L.,  “Vortex  Breakdown  on  Slender 
Sharp-Edged  Wings,”  Journal  of  Aircraft,  Vol.  8,  No.  2,  1971,  pp.  156- 
161. 

“Erickson,  G.  E.,  “Water-Tunnel  Studies  of  Leading-Edge  Vortices,” 
Journal  of  Aircraft,  Vol.  19,  No.  3,  1982,  pp.  442-448. 

*^Hummel,  D.,  “Untersuchungen  uber  das  Aufjplatzen  der  Wirbel  an 
schlanken  Deltaflugeln,”  Zeitschrift  fuer  Flugwissenshaft,  Vol.  13,  1965, 

pp.  158-168.  ... 

^^Khorrami,  M.,  “On  the  Viscous  Modes  of  Instability  of  a  Trailing  Lme 
Vortex,”  Journal  of  Fluid  Mechanics,  Vol.  225, 1991,  pp.  197-212. 

^'^Rockwell,  D.,  Magness,  C.,  Robinson,  O.,  Towfighi,  J.,  Akin,  O.,  Gu, 
W.,  and  Corcoran,  T.,  “Instantaneous  Structure  of  Unsteady  Separated 
Flows  via  Particle  Image  Velocimetry,”  Fluid  Mechanics  Labs.,  Lehigh 
Univ.,  Rept.  PI-1,  Bethlehem,  PA,  Feb.  1992. 

^^Gursul,  I.,  and  Ho,  C.  M.,  “Vortex  Breakdown  over  Delta  Wings  in 
Unsteady  fW  Stream,”  AIAA  Paper  93-0555,  Jan.  1993. 

'^Solignac,  J.  L.,  and  Leuchter,  O.,  “Etudes  experimentales  d’ecoule- 
ments  tourbillonaires  sounds  a  des  effets  de  gradient  de  pression  adverse,” 
Aerodynamics  of  Vortical  Type  Flows  in  Three  Dimensions,  AGARD  CP- 
342,  July  1983. 

'*^Smith,  J.  H.  B.,  “Improved  Calculations  of  Leading-Edge  Separation 
from  Slender,  Thin,  Delta  Wings,”  Proceedings  of  the  Royal  Society,  A, 
Vol.  306, 1968,  pp.  67-90. 

^^Payne,  F.  M.,  Ng,  T.  T.,  Nelson,  R.  C.,  and  Schiff,  L.  B.,  “Visualiza¬ 
tion  and  Wake  Surveys  of  Vortical  Flow  over  a  Delta  Wing,”  AIAA  Jour¬ 
nal,  Vol.  26,  No.  2,  1988,  pp.  137-143. 

^^Schmucker,  A.,  and  Gersten,  K.,  “Vortex  Breakdown  and  its  Control 
on  Delta  Wings,”  Fluid  Dynamics  Research,  Vol.  3, 1988,  pp.  268-272. 

2°Mager,  A.,  “Dissipation  and  Breakdown  of  a  Wing-Tip  Vortex,”  Jour¬ 
nal  of  Fluid  Mechanics,  Vol.  55,  1972,  pp.  609-628. 

^‘Keller,  J.  J.,  Egli,  W.,  and  Althaus,  R.,  “Vortex  Breakdown  as  a  Fun¬ 
damental  Element  of  Vortex  Dynamics,”  Fluid  Dynamics  Research,  Vol. 
3,  1988,  pp.  31-42. 


Appendix  B 


On  fluctuations  of  vortex  breakdown  location 
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Flow  visualization  and  velocity  measurements  are  performed  in  order  to  investigate  the  unsteady 
nature  of  vortex  breakdown  location  over  a  delta  wing.  The  results  indicate  that  the  fluctuations  of 
vortex  breakdown  location  occur  at  much  lower  frequencies  than  the  frequency  of  the  hydrodynamic 
instability  of  the  flow  in  the  wake  of  the  vortex  breakdown.  It  is  suggested  that  the  helical  mode 
instability  does  not  influence  the  unsteady  nature  of  breakdown  location.  ©  1995  American 
Institute  of  Physics, 


It  has  been  observed  in  several  experiments  that  vortex 
breakdown  location  is  not  steady  and  exhibits  fluctuations  in 
the  streamwise  direction.  Garg  and  Leibovich^  reported 
fluctuations  of  the  order  of  the  test-section  radius  in  their 
vortex-tube  experiments.  The  magnitude  of  the  fluctuations 
corresponds  approximately  to  five  times  the  core  diameter 
(after  breakdown)  in  their  experiments.  The  experiments  for 
delta  wings  showed  fluctuations  of  breakdown  location  up  to 
10%  of  the  chord  length  of  the  delta  wing.^  The  purpose  of 
this  Brief  Communication  is  to  report  the  results  of  the  ex¬ 
periments  on  the  fluctuations  of  vortex  breakdown  location. 

One  possible  source  of  the  observed  fluctuations  is  the 
hydrodynamic  instability  of  the  “breakdown  wake”  (i.e.,  the 
flow  in  the  wake  of  the  vortex  breakdown)  where  wake-like 
axial  velocity  profiles  are  observed.  Garg  and  Leibovich^ 
carried  out  velocity  measurements  in  the  breakdown  wakes 
and  observed  coherent  oscillations.  They  suggested  that  these 
oscillations  correspond  to  the  first  helical  mode  instability  of 
the  time-averaged  velocity  profiles.  Periodic  oscillations  due 
to  the  helical  mode  instability  were  observed  in  a  variety  of 
swirling  flows  after  breakdown  occurred.  For  delta  wings, 
the  existence  of  the  first  helical  mode  instability  was  shown 
with  the  help  of  the  two-point  pressure  measurements  on  the 
suction  surface  of  the  delta  wings.^  Examples  of  pressure 
fluctuations^  at  a  fixed  point  underneath  the  vortex  and  near 
the  trailing  edge  are  shown  in  Fig.  1.  At  small  angle  of  attack 
(a=15®),  the  vortex  breakdown  is  not  over  the  wing.  With 
the  increasing  angle  of  attack,  the  vortex  breakdown  moves 
toward  the  apex  of  the  wing.  At  angle  of  attack  a=34®,  the 
breakdown  location  is  approximately  60%  of  the  chord 
length  from  the  wing  apex.  The  length  of  the  time  record  is 
approximately  lOc/Uo,,  where  c  is  the  chord  length  and  I7oo 
is  the  free  stream  velocity.  As  noted  by  Garg  and  Leibovich,^ 
although  the  existence  of  the  wake  instability  has  been  well 
documented,  its  role  in  vortex  breakdown  phenomena  re¬ 
mains  unclear.  The  quasiperiodic  nature  of  the  flow  in  the 
wake  of  breakdown  may  interact  with  the  breakdown  process 
through  upstream  influence.  It  has  not  been  possible  to  de¬ 
termine,  from  the  previous  experiments,  the  role  of  the  heli¬ 
cal  mode  instability  in  the  unsteady  nature  of  breakdown 
location.  The  main  objective  of  this  study  is  to  investigate 
the  relationship  between  the  helical  mode  instability  and  the 
unsteadiness  of  vortex  breakdown  location. 

For  this  purpose,  flow  visualization  and  laser  Doppler 


velocimetry  (LDV)  measurements  were  carried  out  over  a 
delta  wing.  The  experiments  were  performed  in  a  water 
channel  with  a  cross-sectional  area  of  61  by  61  cm.  The  delta 
wing  has  a  sweep  angle  of  A=70®  and  chord  length  of 
c=268  nun.  The  lee  surface  was  flat,  whereas  the  leading 
edges  were  beveled  at  45°  on  the  windward  side.  The  Rey¬ 
nolds  number  based  on  the  chord  length  was  around  Re 
=50  000  and  the  turbulence  level  was  0.6%.  The  velocity 
measurements  were  made  with  a  single  component  LDV  sys¬ 
tem  and  the  free-stream  velocity  was  Uoo-20  cm/s. 

Figure  2  shows  an  example  of  the  streamwise  velocity 
fluctuations  spectrum  taken  just  outside  of  the  vortex  core  at 
an  angle  of  attack  a—3T.  The  velocity  spectrum  at  different 
cross-stream  locations  was  very  similar.  The  axial  location  of 
the  measurement  point  is  approximately  0.08c  downstream 
of  the  time-averaged  breakdown  location.  The  sharp  peak  in 
the  spectrum  is  due  to  the  helical  mode  instability.  The  di¬ 
mensionless  frequency  of  the  instability  fdU^  is  around 
1.72.  This  number  is  in  agreement  with  the  results  from  the 
pressure  measurements^  on  delta  wings,  which  showed  that 
the  dimensionless  frequency  is  on  the  order  of  unity  for  a 
large  range  of  angle  of  attack  and  sweep  angle. 

Flow  visualization  of  vortex  breakdown  was  performed 
by  injecting  fluid  with  food  coloring  dye  near  the  apex  of  the 
model.  The  motion  of  vortex  breakdown  location  was  re¬ 
corded  with  a  30  frames/s  VCR,  a  CCD  camera  and  zoom 
lens.  The  videotape  recording  of  the  motion  was  analyzed 


FIG.  1.  Pressure  fluctuations  near  trailing  edge  (Ar/c=0.93)  for  different 
angles  of  attack.^  Length  of  time  record  is  approximately  lOc/t/*. 
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FIG.  2.  Spectrum  of  streamwise  velocity  fluctuation  in  breakdown  wake, 
vertical  scale  is  linear  and  in  arbitrary  units. 


frame  by  frame  and  the  chordwise  distance  of  breakdown 
location  from  the  apex  was  measured.  The  breakdown  was 
spiral  type,  as  generally  observed  over  delta  wings.  The  lo¬ 
cation  of  breakdown  was  taken  as  the  location  where  the 
streakline  marking  the  core  makes  an  abrupt  kink  to  form  a 
spiral.  The  measurement  uncertainty  for  the  breakdown  loca¬ 
tion  is  0.004c.  The  time  history  of  the  vortex  breakdown 
location  was  observed  for  a  total  length  of  time  record  of 
approximately  lOOc/C/.  (see  Fig.  3).  For  the  experimental 
conditions,  the  time  resolution  is  limited  by  the  frame  speed 
and  conesponds  to  0.025c/CA„.  Therefore,  the  frequencies 
up  to  /c/(7oo=20  can  be  resolved  in  the  frequency  domain. 
For  a=37°,  from  the  time  history  of  breakdown  location,  the 
time-averaged  breakdown  location  and  root-mean-square 
(RMS)  value  of  it  were  calculated  as  i/c -0.297  and 

JCrms/^ 

Figure  4  shows  the  spectrum  of  the  fluctuations  of  vortex 
breakdown  location  for  a-3T  along  with  the  spectra  for 
a=30°  and  26°.  The  frequency  of  the  helical  mode  instabil¬ 
ity  is  also  shown  in  each  plot  with  dashed  lines.  It  is  seen 
that  the  fluctuations  of  breakdown  location  Mcur  at  a  much 
lower  frequency  range  compared  to  the  helical  mode  msta- 
bUity.  It  is  concluded  that  there  is  no  conelation  between  the 
fluctuations  of  breakdown  location  and  helical  mode  instabil¬ 
ity  Indeed,  most  of  the  energy  of  fluctuations  of  breakdown 
location  is  below  /c/f;„=0.2.  For  «=37°,  it  is  also  seen  that 
a  sharp  peak  exists  around  fc/U„^0.12.  The  dominant  fre- 
quency  of  the  oscillations  of  breakdown  location  and  the 
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no.  4.  Spectrum  of  fluctuations  of  breakdown  location;  vertical  scale  is 
linear  and  in  arbitrary  units. 


FIG.  3.  Time  history  of  breakdown  location.  Length  of  time  record  is  ap¬ 
proximately  lOOc/C/uQ. 
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frequency  of  the  helical  mode  instability  vary  with  an^e  of 
attack.  The  fact  that  dominant  frequency  of  the  oscillation  of 
breakdown  location  varies  with  angle  of  attack  suggests  that 
the  osculations  are  not  related  to  the  facUity  dependent  con¬ 
ditions  such  as  the  disturbances  in  the  free-stream  velocity. 

Brief  Communications 


The  most  coherent  oscillations  occur  for  the  largest  angle  of 
attack  a=37®,  when  the  breakdown  location  is  closest  to  the 
apex.  For  all  three  angles  of  attack,  the  main  conclusion  is 
the  same:  the  oscillations  of  breakdown  occur  at  a  much 
lower  frequency  than  that  of  the  helical  mode  instability  in 
breakdown  wake. 

It  should  also  be  noted  that,  the  Kelvin-Helmholtz  in¬ 
stability  of  the  shear  layer  separated  from  the  leading  edge 
has  a  much  higher  frequency  range, /c/C/oo= 10-30. 
Therefore,  it  is  not  related  to  the  fluctuations  of  breakdown 
location  either.  When  compared  with  the  frequency  of  the 
helical  mode  instability,  the  frequency  range  of  the  spectrum 
of  breakdown  location  for  a  stationary  delta  wing  is  much 
closer  to  the  frequency  range  of  typical  aerodynamic  maneu¬ 
vers  (up  to  /c/f/oo=0.03).  The  response  of  breakdown  loca¬ 
tion  and  possible  coupling  between  the  wing  motion  and 
breakdown  location  in  this  frequency  range  is  very  impor¬ 
tant. 
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Active  Control  of  Vortex 
Breakdown  over  a  Delta  Wing 

I.  Gursul,*  S.  Srinivas/  and  G.  Batta* 

University  of  Cincinnati,  Cincinnati,  Ohio  45221-0072 

Introduction 

HE  objective  of  this  study  is  to  accomplish  active  control  of 
vortex  breakdown  over  delta  wings.  The  first  step  in  this  effort 
is  to  identify  a  physical  quantity  that  indicates  the  existence  of  vortex 
breakdown  and  can  be  used  as  a  feedback  signal  for  active  control, 
Pre^ous  studies  suggest  that  pressure  fluctuations  induced  by  the 
helical  mode  instability  of  vortex  breakdown  is  a  good  candidate. 
Measurement  of  pressure  fluctuations  at  a  single  location  on  the 
wing  surface  can  be  sufficient  for  control  purposes.  The  variation 
of  the  amplitude  of  pressure  fluctuations  (or  rms  value  of  pressure) 
with  the  breakdown  location  seems  monotonic.^  Hence,  in  this  study, 
the  rms  value  of  pressure  was  chosen  as  the  control  variable,  and  a 
feedback  control  strategy  was  considered. 

The  second  element  in  active  control  of  vortex  breakdown  is  to 
identify  a  flow  controller  to  influence  the  vortex  breakdown  location. 
Several  methods  were  shown  to  delay  vortex  breakdown.  Blowing 
and  suction  in  the  tangential  direction  along  the  leading  edge,^*^  suc¬ 
tion  applied  around  the  vortex  axis/-^  and  use  of  leading-edge  flaps^ 
are  among  them.  For  most  of  the  techniques  mentioned,  however, 
the  relationship  between  the  control  parameter  and  the  vortex  break¬ 
down  location  is  unknown  or  undesirable  (i.e.,  not  monotonic).  A  de¬ 
sirable  controller  should  have  a  monotonic  relationship  between  the 
control  parameter  and  breakdown  location.  Sweep  angle  has  such  a 
relationship.  Therefore,  variable  sweep  angle  control  was  employed 
in  this  study.  The  relationship  between  the  sweep  angle  and  vortex 
breakdown  location  is  very  well  known  from  static  experiments. 

A  delta  wing  with  variable  sweep  was  fabricated  (see  Fig.  la). 
Measured  rms  value  of  pressure  coefficient  at  a  fixed  point  close  to 
the  trailing  edge  (x/c  —  0.94,  y/c  =  0.27)  is  shown  as  a  function 
of  angle  of  attack  and  sweep  angle  in  Fig.  lb.  The  flat  region  near 
A  =  70  deg  shows  the  pressure  fluctuation  level  in  the  absence  of  the 
vortex  breakdown  over  the  wing  {Cp  =  0.05).  With  increasing  angle 
of  attack  or  decreasing  sweep  angle,  the  vortex  breakdown  moves 
over  the  wing  as  the  rms  pressure  level  increases  and  finally  reaches  a 
saturation.  This  approximate  monotonic  relation  between  the  sweep 
angle  A  and  rms  Cp  suggests  that  a  feedback  control  may  be  feasible. 
The  increase  in  pressure  fluctuations  with  the  decreasing  sweep 
angle  is  due  to  the  increasing  length  of  the  breakdown  region  over  the 
wing,  as  well  as  increasing  circulation  of  the  leading-edge  vortex.^ 
An  important  consideration  for  the  feedback  control  is  the  system 
dynamic  response.  It  is  well  known  that  the  dynamic  response  of 
the  vortex  breakdown  location  in  unsteady  flows  is  characterized  by 
time-lag  effects.  It  was  suggested  that  the  response  of  the  breakdown 
location  is  similar  to  that  of  a  first-order  system.^  Estimated  values 
of  the  time  constant  for  different  types  of  motion  are  summarized 
in  Ref.  8.  Flow- visualization  experiments’  for  variable  sweep  angle 
in  a  water  channel  showed  that  the  normalized  time  constant  is 
tUtsojc  =  2—7.  The  idealization  as  a  first-order  system  suggested 
that  either  proportional  or  integral  control  (or  a  combination  of  both) 
could  be  suitable  for  this  application. 
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A  —  A 

Fig.  la  Schematic  of  variable  sweep  delta  wing. 


Fig.  lb  RMS  value  of  pressure  coefficient  as  a  function  of  angle  of 
attack  and  sweep  angle. 

Experimental  Setup 

Experiments  were  carried  out  in  a  closed-circuit  wind  tunnel  with 
a  cross-sectional  area  of  61  by  61  cm.  Details  of  the  experimental 
setup  and  model  can  be  found  in  Ref.  8.  The  chord  length  of  the 
wing  shown  in  Fig.  la  was  c  =  268  mm  and  the  Reynolds  number 
was  Re  =  190,000.  The  range  of  sweep  angle  was  from  A  =  60  to 
A  =  70  deg.  Two  thin  plates  were  used  to  change  the  sweep  angle. 
The  motion  of  the  plates  was  guided  with  a  cable-pulley-pin  system. 
A  bicycle  brake  cable,  whose  sheath  was  fixed  near  the  trailing  edge 
and  outside  the  wind  tunnel,  was  attached  to  a  drum,  which  was 
driven  by  a  dc  motor  servo  system.  This  flexible  system  allowed  the 
use  of  variable  sweep  even  for  a  pitching  motion  of  the  delta  wing. 
The  unsteady  surface  pressure  was  measured  by  a  high-sensitivity 
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Fig.  2  Block  diagram:  a)  feedback  control  system  showing  signal  conditioning  components  and  b)  with  transfer  functions. 


piezoelectric  transducer  (PCB  model  103 A).  The  measurement  lo¬ 
cation  was  xjc  —  0.94,  y/c  =  0.27,  which  was  chosen  based  on 
previous  measurements.^  The  active  control  experiments  were  con¬ 
ducted  both  for  a  stationary  and  a  pitching  delta  wing.  The  pitching 
mechanism  was  similar  to  the  one  used  by  LeMay  et  al.*^  The  angle  of 
attack  was  varied  between  15  and  23  deg.  A  variable  speed  dc  motor 
and  a  speed  controller  were  used  to  drive  the  pitching  mechanism. 

Since  the  rms  value  of  pressure  was  chosen  as  the  control  variable, 
an  approximation  for  the  rms  pressure  was  necessary  in  a  real-time 
control  loop.  The  raw  pressure  fluctuation  signal  p’ii)  was  rectified 
and  low-pass  filtered.*  The  latter  quantity  is  denoted  by  which 
represents  the  amplitude  of  the  pressure  fluctuations.  The  signal 
Pa  (.0  was  used  as  the  control  variable.  The  overall  objective  of  the 
control  system  was  translated  into  a  control  objective  of  maintaining 
a  particular  value  of  pressure  amplitude  The  implementation 
of  the  controller  can  be  more  easily  understood  by  considering  the 
block  diagram  representation  of  the  system  as  shown  in  Fig.  2a.  The 
plant  in  this  control  system  consists  of  two  parts:  the  mechanical 
position  dynamics  of  the  sweep  angle  mechanism  and  the  vortex 
breakdown  process.  The  pressure  amplitude  /7^(0  was  obtained  by 
rectifying  and  low-pass  filtering,  after  the  signal  from  the  pressure 
transducer  was  passed  through  a  high-pass  filter.  The  rectification 
was  obtained  using  a  full- wave  bridge  rectifier.  The  signal  condi¬ 
tioning  circuits  and  the  feedback  control  system  for  the  sweep  an¬ 
gle  were  implemented  using  operational  amplifiers.  These  circuits 
were  constructed  on  a  breadboard.  The  block  diagram  of  the  feed¬ 
back  control  system  with  the  transfer  functions  is  shown  in  Fig.  2b. 
The  power  amplifier,  dc  motor,  and  sweep  angle  mechanism  were 
modeled  as  a  dc  motor  with  a  directly  coupled  inertia.  Frequency 
response  tests  of  the  sweep  angle  mechanism  indicate  that  this  is 
a  reasonable  approximation.  The  transfer  function  for  this  model 
is  shown  in  Fig.  2b.  The  time  constant  ti  was  estimated  from  the 
step  response  tests  as  Ti  =  0.25  s,  (9,3c/f/oo)«  As  indicated  earlier, 
the  vortex  breakdown  process  is  modeled  as  a  first-order  dynamic 
system  with  time  constant  ti.  The  time  constant  T2  was  estimated 
fern  the  water  channel  experiments*^  as  r2  =  {1-1)c/Uqo, 

Since  the  sweep  angle  positioning  mechanism  provides  an  inte¬ 
gral  relationship  in  the  forward  path,  an  integral  controller  for  the 
overall  control  system  was  relatively  easy  to  implement  In  addi¬ 
tion,  integral  control  has  the  advantage  of  small  steady-state  error. 
The  only  parameter  needed  in  the  controller  design  was  the  gain 
Kp.  Since  specific  values  of  the  plant  parameters  were  not  easily 
obtained,  the  gain  Kp  ^2S  chosen  by  trial  and  error. 

Results 

First,  the  active  control  experiments  were  conducted  for  the  sta¬ 
tionary  delta  wing  at  angle  of  attack  of  =  19  and  23  deg.  The  results 
for  a  =  23  deg  are  shown  in  Fig.  3.  Initially,  the  sweep  angle  was 
set  to  A  =  60  deg  and  the  desired  voltage  uj  (corresponding  to 
a  desired  pressure  amplitude  pAd)  was  chosen  as  1.0  V.  For  this 
value  of  sweep  angle,  the  vortex  breakdown  location  was  over  the 
wing  and  closer  to  the  apex.  In  Fig.  3,  plots  of  pressure  fluctuation 
p\t),  feedback  voltage  Uo(r),  and  sweep  angle  A(r)  are  shown.  At 
around  f  =  5  s,  the  feedback  control  was  turned  on.  Note  that  the 
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Fig.  3  Variation  of  a)  pressure  fluctuation  p\t)y  b)  feedback  voltage 
Va(/),  and  c)  sweep  angle  A(/)  as  a  function  of  time,  a  s  23  deg. 
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Fig.  4  Variation  of  a)  angle  of  attack  a(0,  b)  pressure  fluctuation  c)  feedback  voltage  Va{t),  and  d)  sweep  angle  A(/)  as  a  function  of  time, 
a(/)  =  19  +  4cos(cjf)  (deg),  k  =  uic/Woo  =  0.12. 


sweep  angle  A(r)  increases  and  becomes  nearly  constant  around 
A  =  68.5  deg  after  a  small  overshoot.  This  new  sweep  angle  corre¬ 
sponds  to  a  case  where  the  breakdown  location  is  downstream  of  the 
wing.  The  decrease  in  raw  pressure  fluctuations  and  feedback  volt¬ 
age  accompany  the  variation  of  sweep  angle.  Similar  results  were 
obtained  for  a  =  19  deg.  Experiments  for  different  values  of  the 
desired  voltage  uj,  which  represents  the  desired  pressure  amplitude, 
were  carried  out.  When  this  parameter  was  increased  to  uj  =  1.5  V, 
the  overshoot  of  the  sweep  angle  was  not  observed. 

The  active  control  experiments  for  the  pitching  delta  wing  were 
conducted  in  a  similar  way.  As  the  delta  wing  was  pitched  period¬ 
ically  between  15  and  23  deg,  i.e.,  a{t)  -  19  4-  4cos(4ur)  (deg), 
the  initial  sweep  angle  was  set  to  60  deg.  For  this  range  of  an¬ 
gle  of  attack,  the  vortex  breakdown  was  over  the  wing  during  the 
whole  cycle.  The  plots  of  angle  of  attack  a(r),  pressure  fluctuation 
p'(r),  feedback  voltage  Va(r),  and  sweep  angle  A(r)  are  shown  in 
Hg.  4  for  a  desired  voltage  Vd  =  1.0  Y  and  a  reduced  frequency 
k  =  mcjlU^  =  0.12,  where  m  is  the  radial  frequency  of  the  pitch¬ 
ing  motion.  Again,  the  feedback  control  was  turned  on  suddenly. 
It  is  seen  that  the  sweep  angle  A(r)  increases  and  becomes  nearly 
constant  around  a  new  value  for  which  the  breakdown  location  is 
downstream  of  the  wing.  It  is  also  seen  that  the  magnitude  of  pressure 
fluctuations  decreases  substantially  as  a  result  of  feedback  control. 

Conclusions 

It  is  shown  that  pressure  fluctuations  induced  by  the  helical  mode 
instability  of  vortex  breakdown  can  be  used  to  control  the  vor¬ 
tex  breakdown  location.  Measurement  of  pressure  fluctuations  at 
a  single  location  on  the  delta  wing  is  sufficient  for  this  purpose. 
The  monotonic  variation  of  the  amplitude  of  the  pressure  fluctua¬ 
tions  with  vortex  breakdown  location  makes  the  feedback  control 


possible.  Based  on  previous  experiments  on  the  dynamic  response 
of  vortex  breakdown  location,  the  system  was  idealized  as  a  first- 
order  system,  and  integral  control  was  used.  The  active  control  of 
vortex  breakdown  was  achieved  for  stationary  as  well  as  pitching 
delta  wings. 
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Unsteady  pressure  measurements  were  carried  out  on  a  pitching  delta  wing  in  order  to 
study  the  effects  of  pressure  gradient  on  the  phase  lag  between  the  wing  motion  and  the 
movement  of  breakdown  location.  The  results  indicate  that  the  observed  time  lag  of 
breakdown  location  is  strongly  linked  to  the  external  pressure  gradient  generated  by  the 
wing.  ©  1995  Academic  Press  Limited 


1.  INTRODUCTION 

Several  experimental  studies  of  vortex  breakdown  over  pitching  delta  wings  (see 
Figure  1)  have  shown  that  there  is  a  phase  shift  between  the  motion  of  the  wing  and 
the  movement  of  the  breakdown  location  (Woodgate  1971;  Wolffelt  1987;  LeMay  et  ai 
1990;  Atta  &,  Rockwell  1990).  The  phase  delay  mainly  depends  on  the  reduced 
frequency  (LeMay  et  ai  1990)  and  can  be  as  large  as  180°,  which  means  that  the 
breakdown  location  moves  toward  the  apex  as  the  angle  of  attack  decreases  (Atta  & 
Rockwell  1990).  As  an  example,  the  results  of  flow  visualization  studies  for  sweep 
angle  A  -  70°  are  summarized  in  Figure  2.  Although  the  data  are  for  periodic  pitching 
motion  (see  Figure  1)  with  different  ranges  of  angle  of  attack  (mean  value  and 
amplitude),  there  is  a  consistent  trend  of  increasing  phase  delay  with  increasing 
reduced  frequency.  The  physics  of  this  phase  lag  is  not  well  understood. 

It  is  well  known  that  there  are  two  important  parameters  which  determine  the 
breakdown  location:  swirl  angle  and  external  pressure  gradient  outside  the  vortex  core. 
Hall  (1972)  suggested  that  vortex  breakdown  occurs  as  a  result  of  adverse  pressure 
gradient  along  the  vortex  axis  which  leads  to  a  stagnation  point.  He  showed 
theoretically  that  the  pressure  gradient  along  the  vortex  axis  is  the  sum  of  external 
pressure  gradient  and  swirl  velocity  contribution.  This  is  consistent  with  the  observa* 
tions  that  an  increase  in  swirl  angle  or  in  the  magnitude  of  adverse  external 
pressure  gradient  causes  an  earlier  breakdown  (i.e.,  the  breakdown  location  moves 
upstream).  Hall’s  theory  also  explains  the  observations  that  less  swirl  is  needed  for 
breakdown  if  the  magnitude  of  adverse  external  pressure  gradient  is  increased 
(Sarpkaya  1974). 

For  a  leading-edge  vortex,  the  external  pressure  gradient  is  generated  by  the  wing 
itself,  and  it  is  adverse  due  to  the  presence  of  a  trailing-edge  (Polhamus  1971).  The 
circulation  (or  swirl  angle)  and  adverse  external  pressure  gradient  depend  on 
geometrical  parameters  such  as  angle  of  attack  and  sweep  angle  (Gursul  1992a).  Thus, 
for  a  pitching  wing,  both  the  swirl  angle  and  streamwise  pressure  gradient  vary  during  a 
maneuver. 

It  has  been  suggested  that  delta-wing  vortex  breakdown  in  unsteady  flow  is  driven  by 
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Figure  1.  Definition  diagram  for  the  pitching  delta  wing. 

the  external  pressure  gradient  (Gursul  &  Ho  1994).  Likewise,  the  observed  phase  delay 
of  the  breakdown  location  for  pitching  delta  wings  may  be  related  to  the  variations  of 
the  adverse  pressure  gradient  on  the  wing  surface.  In  order  to  test  this  hypothesis, 
unsteady  pressure  on  a  pitching  delta  wing  was  measured  and  correlated  with  the 
movement  of  the  breakdown  location. 


Figure  2.  Phase  lag  as  a  function  of  reduced  frequency,  A  =  70®:  O,  LeMay  et  al.  (cto  =  34°,  a^=  5°);  •, 
LeMay  et  aL  (oto  =  22*5°,  aj  =  22-5°);  A,  Woodgate  (oq  =  30°,  a j  =  1°);  A,  Woodgate  (ao  =  35°,  -  1°). 
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2.  EXPERIMENTAL  FACILITY 

Experiments  were  carried  out  in  a  closed-circuit  wind  tunnel  with  a  cross-sectional  area 
of  610  mm  by  610  mm.  The  model  delta  wing  (sweep  angle  of  70°)  was  made  of  19  mm 
thick  plexiglas  and  had  a  chord  length  of  c  =  356  mm  (see  Figure  3).  The  lee  surface 
was  flat,  whereas  the  leading-edges  were  beveled  at  45°  on  the  windward  side.  The 
Reynolds  number  was  varied  from  127,000  to  254,000,  although  most  of  the  results  are 
for  Re  =  254,000.  At  the  largest  angle  of  attack,  the  blockage  ratio  was  0-08. 

Unsteady  pressure  measurements  were  made  along  the  centerline  {y/s  =  0)  and 
along  the  ray  y/s  =  0-5.  The  unsteady  surface  pressure  was  measured  by  a  high 
sensitivity  piezoelectric  transducer  (PCB  model  103 A).  The  measurement  uncertainty 
for  the  dimensionless  pressure  coefficient,  Cp,  was  estimated  as  0-005.  The  pressure 
transducer  has  a  built-in  accelerometer  to  cancel  vibration  sensitivity.  It  is  5-6  mm  high, 
9-4  mm  in  diameter  and  has  a  pressure  orifice  measuring  2-54  mm  in  diameter.  The 
sensing  area  of  the  transducer  was  reduced  by  using  a  1  mm  diameter  pinhole  of  2  mm 
depth.  The  amplitude  distortion  and  phase  shift  of  the  pressure  signal  due  to  the  very 
short  pinhole  are  negligible,  since  the  natural  frequency  of  the  system  was  found  to  be 
much  higher  than  the  frequency  range  of  interest.  The  transducer  and  its  leads  were 
housed  in  the  grooves  shown  in  Figure  3.  The  pressure  taps  were  also  used  for  mean 
pressure  measurements.  Additional  pressure  taps  were  located  in  the  spanwise 
direction  at  x/c  =  Q-6  and  x/c  =  0-8.  The  mean  pressure  was  measured  with  a  Setra 
transducer  (model  239).  The  pressure  taps  were  connected  to  the  transducer  with  the 


A 


A 


Cross-section  A-A 


Figure  3.  Top  view  of  the  delta  wing  model  (top),  cross-section  of  the  delta  wing  model  (bottom). 
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Figure  4.  Schematic  of  experimental  set-up. 


enclosed  chamber  (see  Figure  3)  and  plastic  tubing.  All  unused  taps  were  sealed  by 
tape. 

A  schematic  of  the  experimental  set-up  is  shown  in  Figure  4.  The  model  was  pitched 
sinusoidally  about  a  point  located  at  the  one-half  root-chord  position  and  28-5  mm 
below  the  model  suction  surface.  The  pitching  mechanism  was  similar  to  the  one  used 
by  LeMay  et  al.  (1990).  A  variable  speed  DC  motor  and  speed  controller  were  used  to 
drive  the  pitching  mechanism.  The  model  was  pitched  periodically  at  different  reduced 
frequencies  up  to  wc/2t4,  =  0-53.  A  displacement  transducer  was  used  to  monitor  the 
variation  of  the  angle  of  attack.  In  order  to  make  a  comparison  with  the  results  of  flow 
visualization  by  LeMay  et  al.  (1990),  the  angle  of  attack  range  was  chosen  as  29-39°.  In 
addition,  experiments  for  the  range  of  10-15°  (for  which  breakdown  is  absent  over  the 
wing)  were  performed  for  comparison. 

Pressure  signals  were  digitized  and  processed  by  a  laboratory  computer.  The 
dynamic  variation  of  breakdown  location  was  monitored  with  the  help  of  signal  analysis 
of  the  pressure  fluctuations  as  described  in  the  next  section. 

3.  RESULTS 

It  was  shown  in  an  earlier  investigation  that  vortex  breakdown  flowfields  over  delta 
wings  show  helical-mode  instability  which  produces  coherent  pressure  fluctuations  on 
the  suction  surface  (Gursul  1992b).  The  pressure  fluctuations  measured  on  the  steady 
model  (at  x/c  =  0-93,  yls-Q-S)  are  shown  for  different  angles  of  attack  in  Figure  5. 
For  the  smallest  angle  of  attack  a  =  15°,  the  breakdown  location  is  not  over  the  wing. 
For  the  range  of  a  =  29°  to  a  =  39°,  the  breakdown  location  moves  from  xjc  =  0-90  to 
x/c==0-45  (LeMay  et  al.  1990).  With  increasing  angle  of  attack,  the  amplitude  of  the 
pressure  fluctuations  due  to  the  signature  of  breakdown  increases.  This  can  be 
understood  qualitatively  by  modeling  the  vortex  breakdown  as  a  spiral  vortex  filament 
(Jumper  et  al.  1993).  Because  the  circulation  of  the  leading-edge  vortex,  as  well  as  the 
length  of  the  breakdown  region  over  the  wing,  increase  as  the  angle  of  attack  is 
increased,  the  pressure  fluctuations  also  increase  with  angle  of  attack.  Figure  5  shows 
that  the  frequency  content  also  varies  with  the  angle  of  attack.  However,  the  pressure 
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a  ^29^ 


a  =  34“ 


a  =  39“ 


Figure  5.  Pressure  fluctuations  on  steady  model  for  different  angles  of  attack.  Length  of  time  record  is 

250  ms,  xlc  ^0^93,  y/s^O'5. 


signature  of  breakdown  has  a  much  higher  frequency  range  compared  with  the 
frequency  of  a  typical  aerodynamic  maneuver  (Gursul  1994). 

When  the  delta  wing  is  pitched  sinusoidally  (between  29  and  39°),  the  instantaneous 
pressure  signature  of  the  breakdown  varies  as  the  breakdown  location  moves  back  and 
forth  in  the  streamwise  direction  (Figure  6).  It  is  evident  that  the  maximum  fluctuation 
level  due  to  the  breakdown  occurs  not  at  the  maximum  angle  of  attack,  but  later,  at  a 
smaller  angle  of  attack  as  the  wing  pitches  down.  This  is  an  indication  of  phase  delay 
between  the  motion  of  the  wing  and  the  location  of  vortex  breakdown.  In  order  to 
quantify  the  phase  delay,  the  phase-averaging  technique  was  applied  to  the  pressure 
signal,  i.e. 

Cp  =  (Cp>  +  c;, 

where  (Cp)  is  the  phase-averaged  pressure  coefficient.  The  standard  deviation  from  the 
phase-averaged  value, 

was  calculated  as  a  function  of  the  angle  of  attack.  The  variation  of  the  standard 
deviation  shown  in  Figure  7  reveals  hysteresis  loops.  The  loops  become  wider  with 
increasing  reduced  frequency,  which  is  an  indication  of  increasing  phase  delays.  The 
maximum  amplitudes  reached  during  the  maneuver  decrease  in  the  streamwise 
direction  as  the  trailing-edge  is  approached.  The  phase  angles  for  two  streamwise 
locations  are  shown  in  Figures  14  and  15.  These  phase  angles  were  found  by  plotting 
the  variation  of  the  standard  deviation  as  a  function  of  time  (not  shown  here). 
Comparison  of  these  figures  shows  that  the  phase  lag  at  different  streamwise  locations 
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Figure  6.  Pressure  fluctuations  and  time  history  of  angle  of  attack  for  a  periodically  pitching  motion* 
Length  of  time  record  is  838  ms;  xjc  =  0*55,  yjs-  0*5,  reduced  frequency  k  -  0-35. 


Figure  7.  Variation  of  standard  deviation  as  a  function  of  angle  of  attack  at  different  streamwise 
locations  and  yjs-  0-5:  — O — I  k  =  0-19,  — ■ — ,  k  =  0*35,  —A—,  k=Q. 
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is  approximately  the  same,  and  is  in  agreement  with  the  trend  of  the  flow  visualization 
experiments  (Le  May  et  ai  1990). 

In  order  to  investigate  the  possible  relationship  between  the  pressure  field  and  the 
dynamic  character  of  the  breakdown  location,  pressure  measurements  were  carried  out 
along  the  centerline  {yh  =  0)  and  yjs  -  0-5.  For  the  leading-edge  vortex,  the  pressure 
gradient  on  the  suction  surface  of  the  wing  acts  as  the  external  pressure  gradient. 
However,  interpretation  of  the  pressure  along  y/5  =  0-5  is  not  straightforward,  since 
the  measurement  locations  are  approximately  underneath  the  vortex.  Whether  the 
unsteady  variation  of  pressure  is  the  end  result  of  vortex  breakdown  remains  uncertain. 
For  example,  Thompson  et  ai  (1990)  measured  the  pressure  underneath  the  vortex  and 
observed  a  time  lag  with  respect  to  the  wing  motion.  However,  one  cannot  conclude 
whether  the  time  lag  is  because  of  the  movement  of  vortex  breakdown.  On  the  other 
hand,  along  the  centerline  yjs  =  0,  the  effects  of  vortex  on  pressure  are  negligible.  This 
is  verified  by  the  following  experiment.  A  control  cylinder  with  diamater  (i/c  =  0*047 
was  placed  at  a:/c  =  0*92,  z/5  =  0*37  (for  a  =  29°)  which  caused  a  premature  break¬ 
down  (Figure  8).  The  control  cylinder  spanned  across  the  entire  cross-section  of  the 
wind  tunnel.  Although  the  exact  location  of  vortex  breakdown  with  respect  to  the 
cylinder  is  not  known,  its  streamwise  location  is  estimated  to  be  around  xjc  =  0*90 
(LeMay  et  al.  1990)  in  the  absence  of  the  cylinder.  The  breakdown  location  moved 
upstream  after  the  control  cylinder  was  placed.  This  was  monitored  by  the  pressure 
signal  taken  at  the  same  location  as  in  Figure  5.  The  pressure  signals  with  and  without 
the  cylinder  and  the  comparison  with  those  in  Figure  5  confirm  that  the  breakdown 
location  moved  upstream.  The  vertical  location  of  the  cylinder  was  chosen  so  that  the 
largest  amplitude  of  the  pressure  fluctuations  was  obtained.  The  mean  pressure 


Control  cylinder 


Without  cylinder 


With  cylinder 


Figure  8.  Pressure  fluctuations  at  x/c-0'93,  y/s  =  0-5  for  a  =  29°  with  and  without  control  cylinder. 
Dimensionless  diameter  of  the  cylinder  djc  =  0-047;  location  of  the  cylinder  xjc  =  0-92,  zjs  =  0-37.  Dot  on 
the  wing  surface  indicates  the  location  of  pressure  transducer. 
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Figure  9.  Mean  pressure  coefficient  in  spanwise  direction  at  xjc  =  0*80,  a  -  29°;  — O — ,  without  rod; 

— ■ — ,  with  rod. 


coefficient  along  the  spanwise  direction  at  xlc  =  0-80  for  the  cases  with  and  without  the 
cylinder  shows  that  there  is  a  large  reduction  in  the  induced  pressure  underneath  the 
vortex  due  to  the  breakdown,  whereas  the  effects  are  minimum  at  the  centerline 
(Figure  9).  The  measurements  along  the  centerline  (Figure  10)  confirm  that  the 
pressure  is  unaffected  by  the  breakdown,  hence  by  the  vortex  itself. 

The  variation  of  phase-averaged  pressure  along  the  centerline  as  a  function  of  time  is 
shown  in  Figure  11  for  two  different  values  of  the  reduced  frequency.  The  angle  of 
attack  varies  as 

a  =  34°  -  5°  cos  wt. 

The  pressure  along  the  centerline  varies  in  magnitude  as  a  function  of  time.  For  the 
smaller  reduced  frequency  {k  =  0-19),  the  maximum  pressure  (suction)  occurs  when  the 
angle  of  attack  is  closer  to  the  maximum  angle  of  attack;  however,  there  exists  a  time 
delay.  For  k  =  0-35,  this  phase  delay  is  larger.  The  pressure  field  for  a  pitching  wing 
seems  to  be  delayed  in  time  compared  to  the  quasi-steady  case.  Accordingly,  the 
pressure  gradient  has  a  delay  compared  to  the  quasi-steady  case.  This  phase  delay  is 
close  to  that  of  the  breakdown  location  obtained  from  the  pressure  fluctuations  (i.e., 
from  the  variation  of  the  standard  deviation  C'p^  shown  in  Figure  7).  Thus,  there  exists 
a  definite  relationship  between  the  pressure  field  on  the  wing  and  the  location  of  vortex 
breakdown.  Hall  (1972)  showed  that  small  external  pressure  gradients  can  be  amplified 
along  the  core  of  the  vortices,  leading  to  a  stagnation  point.  Thus,  the  large  sensitivity 
of  the  vortex  breakdown  location  to  a  streamwise  pressure  gradient  along  the  exterior 
of  the  vortex  is  very  much  expected.  Recent  computational  studies  (Visbal  1993) 
confirm  that  the  time  lag  of  vortex  breakdown  over  a  pitching  delta  wing  is  linked  to 
the  pressure  gradient  along  the  vortex  axis.  This  pressure  gradient,  which  depends  on 
the  angle  of  attack  and  pitching  motion  (Visbal  1993),  is  strongly  affected  by  the 
external  pressure  gradient  generated  by  the  wing. 
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Chordwise  location,  xlc 

Figure  10.  Mean  pressure  coefficient  in  streamwise  direction  for  yjs  =  0,  a  =  29°:  — O— ,  without  rod; 

— ■ — ,  with  rod. 


The  Fourier  analysis  of  the  phase  averaged  signal  showed  that  most  of  the  energy 
was  concentrated  at  the  fundamental  frequency.  The  phase  angle  between  the  angle  of 
attack  a  and  (-C^)  was  calculated  at  each  measurement  station  and  plotted  in  Figure 
12.  The  phase  delay  increases  with  increasing  reduced  frequency.  At  large  reduced 
frequencies,  there  is  also  considerable  variation  of  the  phase  angle  in  the  streamwise 
direction.  However,  it  should  be  kept  in  mind  that  the  amplitude  of  the  pressure 
variations  becomes  smaller  toward  the  trailing-edge  where  the  largest  phase  angles  are 
observed. 

Pressure  measurements  were  also  made  for  a  smaller  angle-of-attack  range,  for 


Figure  11.  Variation  of  phase-averaged  pressure  along  the  centerline  {yis  =  0)  as  a  function  of  time  for 
a  =  34°  -  5°  cos  wt:  (a)  k  =  0-19;  (b)  k  =  0*35. 
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Figure  12.  Phase  lag  of  phase-averaged  pressure  as  a  function  of  chordv^ise  location  for  a  =  34®  -  5®  cos  cot, 
y/s  =  0.  Re  =  127,000:  — O— ,  k  =  0*19,  -d—,  k  =  0-36,  -A--,  k  =  0*53. 

which  the  breakdown  is  absent  (this  was  checked  by  monitoring  the  pressure 
fluctuations): 

a  =  12*5®  “  2*5®  cos  cut. 

Similar  variations  of  phase  lag  shown  in  Figure  13  indicate  that  the  source  of  the  phase 
delay  is  not  the  breakdown  process  itself.  Similar  phase  lags  for  the  normal  force 
coefficient  were  predicted  with  a  panel  method  in  the  absence  of  vortex  breakdown 
over  a  delta  wing  (Ashley  et  al  1991).  Thus,  there  is  enough  evidence  that  the  external 
pressure  gradient  is  responsible  for  the  observed  phase  lag  of  breakdown  location. 

For  the  case  with  breakdown,  i.e.  a  =  34°  -  5°  cos  u)t,  the  variation  of  phase  lag 
between  the  wing  motion  and  {-Cp)  at  two  streamwise  locations  is  shown  as  a  function 
of  reduced  frequency  in  Figures  14  and  15  for  different  values  of  Reynolds  number. 
The  results  show  that  the  measured  phase  lags  are  not  sensitive  to  the  Reynolds 
number.  Also  shown  are  the  phase  lags  obtained  from  the  pressure  fluctuations  (i.e., 
the  variation  of  the  standard  deviation  Cp^^  shown  in  Figure  7)  and  determined  by  flow 
visualization  for  a  similar  model  (LeMay  et  al.  1990). 

4.  CONCLUSIONS 

In  order  to  study  the  effect  of  pressure  gradient  on  vortex  breakdown,  unsteady 
pressure  measurements  on  a  pitching  delta  wing  were  carried  out.  It  was  shown  that  the 
pressure  fluctuations  induced  by  the  helical  mode  instability  of  vortex  breakdown  can 
be  used  to  quantify  the  phase  lag  between  the  wing  motion  and  breakdown  location. 
Variation  of  the  standard  deviation  of  pressure  fluctuations  from  the  phase-averaged 
value  is  shown  to  be  a  useful  tool  to  detect  the  vortex  breakdown. 

With  the  help  of  time-averaged  pressure  measurements  with  and  without  forced 
breakdown,  it  was  shown  that  the  pressure  along  the  centerline  {y  Is  =  0)  (where  the 
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Figure  13.  Phase  lag  of  phase-averaged  pressure  as  a  function  of  chordwise  location  for  a  =  12'5'’  - 
2-5°  cos  a)t,  y/s  =  0,  Re  =  127,000;  — O— ,  k  =  0-19,  — ,  k  =  0-35,  —A—,  k  =  0-52. 

effect  of  the  vortex  is  negligible)  is  representative  of  the  external  pressure  field  for  the 
leading-edge  vortices.  Variation  of  the  phase-averaged  pressure  along  the  wing 
centerline  as  a  function  of  time  shows  that  there  is  a  phase  delay  between  the  wing 
motion  and  the  pressure  field.  This  phase  delay  is  close  to  the  phase  delay  of  the 


Reduced  frequency 

Figure  14.  Phase  lag  as  a  function  of  reduced  frequency  for  x/c  =  0*72,  a  =  34®  -  5®  cos  a>/:  O,  <-C^>, 
Re  =  127,000;  □,  Re  =  190,500;  A,  <-C^),  Re  =  254,000;  A,  Cp^;  ■,  visualization  (LeMay  et  al.  1990). 


582 


L  GURSUL  AND  H.  YANG 


Figure  15.  Phase  lag  as  a  function  of  reduced  frequency  for  x/c  =  0-89,  a  =  34®  -  5®  cos  O, 

Re  =  127,000;  □,  {-Cp),  Re  =  190,500;  A,  <-Cp),  Re  ==  254,000;  A,  ■,  visualization  (LeMay  et  ai  1990). 

movement  of  breakdown  location.  The  phase  delay  of  pressure  increases  with 
increasing  reduced  frequency  and  is  not  sensitive  to  Reynolds  number  variation. 
Similar  phase  lags  for  pressure  were  found  for  a  smaller  angle-of-attack  range  for  which 
the  breakdown  is  absent.  This  confirms  that  the  observed  time  delays  in  breakdown 
location  is  related  to  the  variation  of  the  external  pressure  gradient  generated  by  the 
wing. 

The  phase  lag  of  the  external  flow  field  is  due  to  the  shedding  of  vorticity  from  the 
wing  (Ashley  et  al  1991).  It  should  also  be  noted  that  the  swirl  level,  which  is  the  other 
important  parameter  for  vortex  breakdown,  was  not  measured  in  this  study.  However, 
the  pressure  measurements  suggest  that  the  external  pressure  gradient  plays  a  major 
role. 
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APPENDIX:  NOMENCLATURE 

chord  length 

pressure  coefficient 

phase-averaged  pressure  coefficient 

fluctuating  pressure  coefficient 

standard  deviation  of  pressure  coefficient 

diameter  of  control  cylinder 

reduced  frequency,  a)c/2L/* 

pressure  fluctuations 

Reynolds  number,  LLc/v 

local  semispan 

freestream  velocity 

chordwise  distance  from  wing  apex 

spanwise  distance  from  wing  root 

distance  above  wing  surface 

angle  of  attack 

mean  angle  of  attack 

amplitude  of  pitching  motion 

sweep  angle 

radial  frequency 

phase-averaged  quantity 
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Abstract 

Effect  of  leading-edge  devices  on  vortex 
breakdown  was  investigated.  Stationary  as  well  as 
oscillating  leading-edge  flaps  were  shown  to  alter  the 
structure  of  the  leading-edge  vortices,  and 
consequently  the  location  of  vortex  breakdown.  Also, 
periodic  variations  of  the  sweep  angle  were  used  to 
control  the  breakdown  location  for  a  pitching  delta 
wing. 

Introduction 

Control  of  vortex  breakdown  continues  to  be 
of  vital  importance  because  the  breakdown  may  have  a 
considerable  effect  on  aircraft  performance,  such  as 
the  effect  on  the  time-averaged  lift  force*.  Breakdown 
is  also  important  because  of  the  unsteady  nature  of 
flow  downstream  of  vortex  breakdown^.  The 
unsteadiness  may  affect  the  stability  of  the  aircraft  and 
also  cause  buffeting.  These  concerns  have  stimulated 
research  on  potential  (active  and  passive)  control 
methods  to  delay  vortex  breakdown. 

This  control  problem  is  particularly  important 
for  unsteady  delta  wings  as  well  as  for  unsteady  free 
stream^'^.  The  vortex  breakdown  may  appear  on  the 
wing  even  if  no  breakdown  is  observed  in  the  steady 
case®.  The  time-dependent  nature  of  breakdown 
location  may  affect  unsteady  loading  on  the  wing. 
One  of  the  objectives  of  this  work  is  to  study  the 
possibility  of  controlling  the  breakdo^vn  in  unsteady 
flows. 

It  is  well  known  that  there  are  two  important 
parameters  which  determine  the  breakdown  location; 
swirl  angle  and  external  pressure  gradient  outside 
the  vortex  core.  The  principle  of  any  control  method 
is  to  alter  either  or  both  of  these  parameters.  Blowing 
and  suction  in  the  tangential  direction  along  a  rounded 
leading-edge’’*,  and  suction  applied  around  the  vortex 
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axis®’*®  were  shown  to  delay  breakdown.  The  former 
is  believed  to  be  due  to  a  change  in  the  strength  and 
location  of  the  leading-edge  vortices.  This  is  achieved 
by  affecting  the  location  of  separation  on  the  rounded 
leading-edge.  Direct  control  of  the  separation  point 
modifies  the  vortex  properties  such  that  a  reduction  in 
the  swirl  angle  is  achieved.  The  piupose  of  applying 
suction  arotmd  the  vortex  axis  is  to  reduce  the  local 
adverse  pressure  gradient,  which  is  achieved  by 
accelerating  the  axial  flow  along  the  core. 

Since  all  of  the  vorticity  of  the  leading-edge 
vortices  originates  from  the  separation  point  along  the 
leading-edge,  leading-edge  devices  are  particularly 
attractive  tools  that  can  be  used  to  influence  the 
strength  and  structure  of  these  vortices.  For  example, 
leading-edge  flaps  are  known  to  be  capable  of 
controlling  the  circulation  and  location  of  the  leading- 
edge  vortices**'**.  These  studies  concentrated 
primarily  on  the  effect  of  leading-edge  flaps  at  low 
angle  of  attack  where  vortex  breakdovm  was  not 
observed  over  the  wing.  Recently,  the  effect  of  vortex 
cavity  flaps  at  high  angles  of  attack  has  been 
reported*  “*.  The  purpose  of  this  paper  is  to  explore  the 
effect  of  leading-edge  flaps  on  the  vortex  breakdown 
phenomena.  Stationary  as  well  as  oscillating  flaps 
may  provide  further  insight  into  the  possibility  of 
control  of  vortex  breakdown.  Also  studied  was  vortex 
breakdown  over  a  delta  wing  with  variable  sweep 
angle,  since  this  method  allows  a  direct  control  of  the 
strength  and  structure  of  the  leading-edge  vortices. 
Again,  periodic  variations  of  sweep  angle  were  studied 
with  the  control  of  vortex  breakdown  in  mind. 
Application  of  this  technique  in  unsteady  flows  was 
demonstrated  for  a  pitching  delta  wing. 

Experimental  Facility 

Experiments  were  carried  out  in  a  water 
charmel  with  a  cross-sectional  area  of  61  cm  by  61  cm. 
The  model  delta  wing  with  leading-edge  flaps  and  the 
variable  sweep  delta  wing  are  shown  in  Figures  1  and 
2,  respectively.  For  both  models,  the  base  model  has  a 
sweep  angle  of  A=70*.  The  chord  lengths  are  c=254 
mm  and  268  nun  for  the  model  with  flaps  and  the 
variable  sweep  delta  wing,  respectively.  The  Reynolds 
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breakdown  appears  on  the  wing  when  the  swirl  level 
reaches  a  certain  value.  For  a=20°,  the  swirl  level 
does  not  change  much  with  the  flap  angle  5  when  the 
breakdown  is  over  the  wing.  On  the  other  hand,  for 
a=25°,  the  swirl  level  increases  rapidly  with  the  flap 
angle  (when  the  breakdown  is  over  the  wing), 
explaining  the  sensitivity  of  breakdown  location  with 
respect  to  the  variations  in  the  flap  angle. 

In  addition  to  the  changes  in  the  swirl  level 
(or  strength  of  the  vortex),  the  location  of  the  vortex 
core  (in  the  cross  plane)  exhibits  variations  with  the 
flap  angle.  The  location  of  the  vortex  core  was 
obtained  from  the  velocity  measurements  (not  shown 
here),  which  suggested  that  the  existence  of  the  flap 
may  affect  the  pressure  gradient  outside  of  the  vortex 
core.  For  example,  for  a=20°,  there  is  a  gradual 
variation  of  breakdown  location  with  the  flap  angle, 
although  the  swirl  level  does  not  change  much  in  the 
same  range.  This  implies  that  the  external  pressure 
gradient  for  the  vortex  core  is  affected  by  the  changes 
in  the  geometry  of  the  wing. 

E.xperiments  with  oscillating  flaps  were  also 
conducted.  The  variation  of  vortex  breakdown 
location  for  harmonic  variations  of  flap  angle  between 
60®  and  180°  is  shown  in  Figure  7  together  with  that 
of  the  steady  flap  angles.  The  reduced  frequency  is 
k=a)c/2Ua,=0.4,  where  cn  is  the  radial  frequency  and  U« 
is  the  free  stream  velocity.  The  variation  of 
breakdown  location  reveals  a  hysteresis  loop.  This  is 
an  indication  of  a  time  lag  response,  which  is  similar 
to  the  response  of  a  pitching  delta  wing.  Clearly, 
development  of  the  detailed  flow  structure  is  needed  in 
order  to  understand  the  response  of  vortex  breakdown. 
For  this  purpose,  experiments  to  obtain  the  phase- 
averaged  flow  field  for  harmonic  variations  of  flap 
angle  are  currently  underway. 

Vortex  Breakdown  over  a  Delta  Wing  with 
Variable  Sweep  Angle 

In  static  experiments,  the  location  of  vortex 
breakdown  moves  toward  the  apex  with  decreasing 
sweep  angle  (or  aspect  ratio),  as  Imown  from  a  number 
of  e.xperiments  with  different  aspect  ratio  wings. 
Therefore,  in  the  steady  case,  the  location  of 
breakdown  shows  a  monotonic  variation  with  the 
sweep  angle,  as  opposed  to  the  variation  with  the  flap 
angle.  This  feature  may  be  advantageous  for  active 
control  purposes'*. 

The  variation  of  breakdown  location  was 
studied  for  harmonic  variations  of  s\veep  angle,  i.e. 

A=Ao+Ai  sin  ojt 

The  variations  of  the  ensemble-averaged  breakdown 
location  are  shown  in  Figure  8  for  Ao=65°  and 


Ai=2.5°  for  different  values  of  the  reduced  frequency, 
which  show  hysteresis  loops.  The  loops  become 
narrower  with  increasing  reduced  frequency. 
Actually,  at  very  large  reduced  frequency  (k=2.0  and 
k=4.0  (not  shown  here)),  the  ensemble-averaged 
breakdown  location  is  almost  constant.  The  phase  lag 
of  vortex  breakdown  location  with  respect  to  that  of 
the  quasi-stea^  case  was  found  by  plotting  the  time 
history  of  breakdown  location.  This  is  shown  in 
Figure  9.  It  seems  that  the  phase  lag  reaches  a 
maximum  at  k=0.4  and  then  decreases.  This  time-lag 
behavior  seems  to  be  an  inherent  response  of 
breakdown  location  in  unsteady  flows,  regardless  of 
the  type  of  motion'*’'*. 

In  order  to  understand  the  phase  lag  of  vortex 
breakdown,  phase-averaged  velocity  measurements 
were  carried  out  in  a  cross  plane  (at  x/c=0.5)  upstream 
of  breakdown  location  for  a  reduced  frequency  of 
k=0.4.  It  should  be  noted  that  the  flow  field  of  the 
leading-edge  vortex  develops  very  gradually  before 
vortex  breakdown.  Indeed,  the  velocity  changes 
drastically  only  in  a  small  neighborhood  of  the 
breakdown  location.  It  is  also  known  that  the  swirl 
level  of  the  leading  edge  vortex  does  not  vary  much  in 
the  streamwise  direction*.  Therefore,  these 
measurements  give  a  good  indication  of  the  upstream 
conditions  before  breakdown.  Contours  of  constant 
phase-averaged  axial  velocity  are  shown  at  different 
sweep  angles  for  Ao=66®  and  Ai=3°  in  Figure  10. 
Also  shown  are  the  mean  velocity  contours  in  the 
steady  case  for  the  minimum  sweep  angle  A=63°  and 
the  maximum  sweep  angle  A=69°.  As  the  sweep 
angle  varies  dynamically,  the  location  of  the  vortex 
core  (defined  as  the  location  at  which  the  phase- 
averaged  axial  velocity  is  maximum)  varies  mostly  in 
the  spanwise  direction.  Note  that  the  axial  velocity 
contoius  are  very  similar  for  the  dynamic  and  static 
cases.  Also  the  two  flow  fields  for  the  mean  sweep 
angle  Ao=66®  (for  increasing  and  decreasing  sweep 
angle)  are  not  very  different  from  each  other.  All 
these  suggest  that  the  structure  of  leading  edge  vortex 
is  not  affected  much  by  the  unsteadiness  and  is 
approximately  quasi-steady.  Nevertheless,  there  is  a 
sniall  phase  lag  in  the  development  of  the  vortical 
flow.  This  is  evident  from  the  contours  of  constant 
standard  deviation  from  the  phase-averaged  velocity 
shown  in  Figure  11.  It  should  be  noted  that  the  shape 
of  the  contours  of  the  standard  deviation  (or  rms 
velocity  fluctuations  in  the  steady  case)  is  related  to 
the  main  vortex  as  well  as  the  shear  layer  and 
secondary  separation  region. 

The  variations  of  the  main  parameters  are 
summarized  in  Figure  12.  The  variation  of  sweep 


Conclusions 

Control  of  vortex  breakdown  with  leading- 
edge  devices  was  considered.  The  main  idea  is  to 
control  the  strength  and  location  of  leading-edge 
vortices,  which  in  turn  determine  the  swirl  level  (and 
possibly  the  external  pressure  gradient)  for  the 
swirling  flow. 

The  effect  of  leading-edge  flaps  on  vortex 
breakdown  was  investigated.  It  was  shown  that  the 
effect  of  flaps  and  sensitivity  of  breakdown  location 
strongly  depends  on  angle  of  attack.  Velocity 
measurements  showed  that,  with  the  varying  flap 
angle,  large  changes  take  place  in  the  vortex  core 
diameter  and  location,  in  axial  and  swirl  velocity 
profiles,  and  in  the  secondary  vortex  structure.  In 
general,  with  the  increasing  flap  angle,  larger  swirl 
velocities  are  observed.  The  variation  of  maximiim 
swirl  velocity  seems  to  be  correlated  with  the 
variations  of  breakdown  location  as  the  flap  angle  is 
varied. 

E,xperiments  with  oscillating  flaps  showed 
that  there  e.xists  a  hysteresis  loop  in  the  variation  of 
breakdown  location.  This  is  an  indication  of  a  time- 
lag  response,  which  is  similar  to  the  response  of  a 
pitching  delta  wing. 

Harmonic  variations  of  sweep  angle  were 
studied  for  a  delta  wing  with  variable  sweep.  The 
variation  of  vortex  breakdown  location  showed 
hysteresis  loops  and  phase  lag  which  depend  on  the 
reduced  frequency.  However,  measurements  of  the 
phase-averaged  axial  velocity  in  a  cross  flow  plane 
upstream  of  breakdown  location  showed  no  major 
unsteady  effects  in  the  development  of  the  flow  field. 
The  larger  time  lag  of  vortex  breakdown  location 
caimot  be  explained  with  the  variations  in  the 
upstream  conditions  since  the  location  and  strength  of 
the  vortex  exhibit  smaller  phase  lags.  This  suggests 
that  the  external  pressure  gradient  plays  a  major  role 
in  the  dynamic  response  of  breakdown  location. 

Use  of  variable  sweep  for  breakdown  control 
was  demonstrated  for  a  pitching  delta  wing. 
Oscillations  of  sweep  angle  with  the  same  frequency  as 
pitching,  but  with  a  phase  angle  were  considered.  It 
was  shown  that,  at  an  optimum  phase  angle,  the 
perturbations  of  sweep  angle  decrease  the  amplitude  of 
the  variations  of  breakdown  location.  In  some  cases,  it 
was  also  observed  that  the  average  breakdown  location 
moved  downstream. 
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Figure  3:  Schematic  of  experimental  setup. 
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Figure  4:  Variation  of  breakdown  location  as  a  function  of  flap  angle  for  several  values  of  angle  of  attack. 
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Figure  5:  Constant  contours  of  normalized  mean  axial  velocity  for  5=0“  (top,  a=20‘’),  5=60“  (middle 
a=25“),  5=100“  (bottom,  a=20“). 


Figure  6:  Variation  of  maximum  swirl  velocity  as  a 
0  function  of  flap  angle  for  a=20°  and  ct=25°.  The 

arrows  indicate  the  flap  angle  at  which  breakdown  is 
at  the  trailing-edge  of  the  wing. 


Figure  7:  Variation  of  vortex  breakdown  location  for 
harmonic  variations  of  flap  angle  (k=0.4),  and  for 
static  flaps  ((x?=30®). 


Figure  9:  Variation  of  phase  lag  as  a  function  of 
reduced  frequency. 
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Figure  8:  Variation  of  breakdown  location  for 
harmonic  variations  of  sweep  angle  and  for  static  case. 
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Figure  13:  Variation  of  breakdown  location  for 
pitching  motion  and  for  static  case. 


Figure  15:  Maximum  and  minimum  locations  of 
breakdown  as  a  function  of  phase  angle.  The  dashed 
line  is  for  pitching  motion  only  for  the  average  sweep 
angle. 


Figure  16:  Maximum  and  minimum  locations  of 
breakdown  as  a  function  of  phase  angle.  The  dashed 
line  is  for  pitching  motion  only  for  the  average  sweep 
angle. 


Figure  14:  Variation  of  breakdown  location  for 
combined  motion  of  pitching  and  variable  sweep  for 
different  values  of  phase  angle. 
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